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1.0 IMTRODUCTION AND PURPOSE 

This document presents a comprehensive survey of past and current research 
efforts in the area of geometrically nonlinear finite elements. The survey 
is intended to serve as a guide in the choice of nonlinear elements for 
specific problems, and as background to provide directions for new element 
developments. The elements are presented in a handbook format and are 
separated by type as beams, plates (or shallow shells)., shells, and other 
elements. Within a given type, the elements are identified by the assumed 
displacement shapes and the forms of the nonlinear strain equatiohs. Solu- 
tion procedures are not discussed except when a particular element formula- 
tion poses special problems or capabilities in this regard. 

The main goal of the format is to provide quick access to a wide variety of 
element types, in a consistent presentation format, to facilitate comparison 
and evaluation of different elements with regard to features, probable 
accuracy, and complexity. References are provided to aid in further detailed 
studies which readers may wish to undertake. 

The survey is based primarily on published literature. In addition, however, 
in order to have an up-to-date list, researchers have been contacted by 
letter and requested to prepare inputs to this document based on their most 
current efforts. 
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2.0 kOPE 

The available literature on nonl infear finite elements is quite large. The 
types of elements, the basis of derivations and solutions, and the types of 
nonlinearities considered are all extensively documented in the literature. 

Some limitations have been imposed on this document in order to obtain a 
useful yet comprehensive compilation of nonlinear elements. 

First, emphasis has been placed on geometric nonlinearities as opposed to 
material nonlinearities. Geometric nonlinearities involve basic assumptions 
concerning strain-displacement equations and the development of new elemental 
matrices. They also strongly Influence the relative adequacy of element dis- 
placement functions. Thus, the geometric nonlinearities of the formulations 
are intrinsic to the elements, while extensions to nonlinear material char- 
acterization can be incorporated In most elements by employing numerical 
Integration schemes in calculations of element properties, and by using 
highly competent solution procedures. For these reasons, the research for 
this document has been directed toward geometric non 1 i near i t les and material 
non 1 inear i t ies have only been included as an adjunct to this goal. 

Secondly, only finite elements based on the displacement method have been 
considered. Neither the force or equilibrium elements, nor the various 
hybrid elements have been included. This restriction Is motivated by the 
fact the hybrid and force method elements employ a considerably different 
approach than displacement elements, and are difficult to classify In the 
latter group . 

Third, solution methods have not been emphasized in the document. This subject 
has been reviewed rather extensively in the published literature. Occasionally, 
the derivation of elemental matrices Is tied closely to the solution procedure 
used to solve the structural problem. In these cases the usefulness or accuracy 
o‘' the element may not be immediately evident by simply considering elemental 
matrices. However, in general, a variety of solution procedures can be used 
with any given element formulation. Thus, the omission of this subject herein 
does not compromise the usefulness of the element descriptions given, in 
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general, a variety of solution procedures can be used with any given 
element formulation. Thus, the emission of this subject herein does not 
compromise the usefulness of the element descriptions given. In some cases 
where solution procedures are specialized for a certain element, they are 
discussed briefly. 

Qualitative comparison of the various elements has not been made. Those 
advantages and disadvantages and relative accuracy data that have been 
noted in the literature are repeated where the data appear particularly 
important. 

Finally, the element descriptions in the handbook include elements for 
which complete formulations and examples of problem solutions have been 
published. There are a number of very general finite element formulations 
which have not been specialized to apply to particular elements, but which 
could be used to generate a number of related elements. Such general 
formujations are described in References 10 ^ through 107 - The last refer- 
ence is particularly complete in this regard. 
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3.0 USER'S GUIDE 


-The handbook is composed of four parts based on element type. These 
basic types are: 


B-1,3, ^ 

P-1,2,- 

S-1,2, 

C-1,2,— — 


Beam or line elements 

Plate or shallow shell elements 

Shdll Elements 

Solid or three-dimensional elements 


Within a given element type, the elements are catalogued by displacement 
function. Each new displacement function is given a new element I.D. (i.e., 
P-2, P-3, P-4). When several authors have used the same displacement func- 
tion for a given element but differences in scope, features, and capabilities 
exist between the resulting elements, then the element identification retains 
the basic I.D., but carries a lower case alphabetic delineator (i.e., B-2a, 
B-2b, B-2c) . The handbook description given for each. element includes the 
I.D. number and a brief definition of element type, followed by the displace- 
ment function equations and concise description of the element in narrative 
form. This is followed by the reference where the element was derived, 
variations from the basic element, and any known advantages or disadvantages 
of the element. The strain displacement equations are given to indicate how 
the geometric nonlinearities are introduced into the strains, and what non- 
linearities are included in the formulation. This is followed by a discussion 
of the element in which user-orientated information concerning material 
description, coordinate systems, and solution procedures are presented. 

No attempt has been made to list chronologically the authors of various 
finite element concepts , In many cases several authors have published slight 
variations of a single element, and it appears preferable to deal with the 
elements themselves rather than the chronology of development. The order 
of the references given in this document is thus quite random. 

In many cases the plate elements are applicable to shallow shells. The class- 
ification of such elements overlaps to some extent. It is therefore suggested 
that the reader refer to both the section on plates and the section on shells 
for each application. 


4.0 NOMENCLATURE 


In forming a reference book of this kind, based solely on the published 
works of many authors, a great variety of nomenclatures were encountered. 
Early in this project the decision was made not to convert the variety of 

nomenclatures to a common nomenclature (a formidable task) but rather to use 

the reference authors own nomenclature in the writeup for that reference. 

The reason being that the purpose of this book is not to present programmable 

matrices but rather to gather together the published nonlinear elements and 
provide sufficient inforraatiCn for a reader to direct his attention to a 
particular reference of interest. The reference material is obviously the 
best source for a detailed study of the elemental matrices, elemental 
characteristic, solution procedure, and applications. 


Much of the nomenclature is common in the various references, such as: 


E 

a 

E 

IK] 

IK°1 

IK'] 

{ 6 } 

u,v,w 

x,y,z 

DOF 


- Young’s Modulus 

- stress 

- strain 

- stiffness matrix 

- linear stiffness matrix 

- K - one geometric stiffness matrix 
column matrix of nodal displacements 

- displacements 

rectangular cartesian coordinates 

- degrees of freedom 


When an author uses unique symbols or nonstandard nomenclature these terms 
are defined arid in addit-’on, one or more figures are generally used to 
illustrate the authors piirticular notation. Subscripts are often used to 
denote derivatives. Their use is obvious in most applications. 
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5.0 NONLINEAR STRAINS, COORDINATE SYSTEMS, AND SOLUTION PROCEDURES 

This section briefly discusses the sequence of calculations , particularly 
those associated with stress and stiffness matrix calculation, which is 
involved in nonlinear finite element problem solutions. The goal of this 
discussion is to offer some perspective concerning the several ways of 
handling the nonlinear strain displacement equations. The subjects dis- 
cussed here are a part of the overall subject of solution procedures for 
nonlinear finite element problems. Several excellent references are avail- 
able which discuss this subject in more detail and from different viewpoints 
from those used here. References 10 A through 107 are particularly helpful 
in this regard. 

Global Coordinate System and Nonlinear Strain Equations 

Finite element derivations nearly always adhere to the basic Lagrangian 
formulation of mechanics In this formulation the displacements and forces 
are referred to a fixed, global coordinate system. The strains and stresses 
are "convected" with the material as the deformation progresses and strains 
are computed by referring distortions to the initial lengths and orientations 
of the material elements. For a cartesian global coordinate system, the 
strain of the material "fiber" initially aligned with the x axis is given by 
the familiar formula, which occurs frequently in the element descriptions 
in this report. 



The incremental strain is given by, 


£i£v = 





2> V 

bv: 






"bx 'hx 






e is actually a component of the so-called second Piola-Kirchhof f strain 
tensor, rather than a true physical Strain. However, for small strains, 
for cartesian coordinate systems, it can be taken as the physical 
strain with sufficient accuracy. 

Figure 1 illustrates for a very simple case two facts concerning the 

application of the above equations which are particularly important for 

finite element calculations. In the figure, the w and u displacements are 

both large and ( ~ )^and ( ~ are of roughly equal importance in the 

3u 

calculation of the strain even though the strain is small. Note that 

is not a good approximation to the strain. In the early stages of the 

movement from the initial to the final positions of the bar element, 

3u 

however, is small and is a good approximation to the strain. Hence, 

early in the deformation, i.e., while the bar is still close to parallel 

9u 2 

to the X axis, the term ( ) can be omitted in the strain equation, 

9w 2 

though ( — ) must be retained. 

Figure 1 also illustrates the fact that physical meanings of the u and w 
displacements change as the bar rotates away from the position parallel to 
the X axis. Early in the deformation, u is essentially an axial displace- 
ment, and w is a lateral, or bending, type of displacement. This is signif- 
icant for finite element work, because usually the element displacement 
functions used for the axial and bending types of displacements differ from 
one another, with the bending displacement functions being considerably more 
competent than the axial functions. In the latter stages of the deformation 
illustrated by the figure, the u displacement has the nature of a bending 
displacement. If simple element functions were used for u, they would 

thus become Incompetent for large deformations. 

If the global coordinate system is curvilinear, as it might be for a shell 
problem, the nonlinear strain equations become more complicated. The 
essential features of the above discussion remain unchanged, however. 
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A final point needs to be made concerning the nonlinear strain equation. 

Figure 1 shows both solid-line and dashed-line positions of the bar, cor- 
responding to non-bending and bending types of problems, respectively. 

The approximation is frequently madej for bending problems, that the non- 
linear terms in the strain equation only include the straight-line element 
shape, i.e., ignore the bending curvature. Thus, the use of the term 
( 1"^ in an element formulation is not completely definitive; the w may 
or may not include the element "bowing" during the actual calculation of the 
strain. 

Stiffness Matrices, Incremental Solutions 

If the nonlinear strain equations are substituted into an energy or virtual 
work integral to derive the element incremental equilibrium equations, a 
number of different linear and nonlinear terms are obtained. Defining Aq, 

K, and AP, respectively, to be element incremental displacements, stiffness 
matrices, and incremental loads, the linearized matrix incremental quilibrium 
equation has the form 

(K + K„)Aq = L? 
o (j 

The Ko matrix contains the stiffness of the element against very small de- 
formations caused by the incremental displacements. is not a constant 

matrix, since it includes the effects of accumulated deformations prior to 
the increment in question. can be thought of as containing two parts, 
the first representing linear behavior associated with very small displace- 
ments from the initial, undeformed state; and the second accounting for the 
changes in stiffness due to the accumulated rotations. Often only the first 
part is used in an element formulation, resulting in a constant matrix. 

K contains stresses, and provides an "effective" stiffness effect associated 
G 

with the changes of direction of the stresses in the element as the rotations 
due to Aq take place. The matrix, often called the geometric, differential, 
or initial stress stiffness matrix, arises from the nonlinear terms in the 
strain equation. The non-constant part of arises from these same terms, 

and the linearization to obtain a constant matrix therefore 
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results from dropping the same nonlinear terms which gave the matrix. 

The above linearized Incremental equilibrium equation has Itself resulted 
from dropping second order terms In Aq, these terms also arising from non- 
linear terms In the strain equations. 

It Is seen that an element's use of nonlinear strain equations, such as 
those given in the element descriptions In this handbook, does not necessarily 
mean that the nonlinear element In question will In calculations demonstrate 
fully nonlinear behavior. The nonlinear character of the element depends 
on how the nonlinear strains equations are used In the element derivation 
and In problem solutions. Several options are used In this regard, partly 
associated with solution procedures, and partly with the element equations 
themselves . 

The NASTRAN program uses a simple but often effective procedure In which 
no nonlinear strain development Is Included at all, but the effect Itself 
is retained. In this case the strain and stress development Is based on the 
linearized stra In-d Isplacemqnt equations, the.K^ matrix Is the constant 
initial state matrix, and through the matrix the nonlinear effects of the 
rotations on the equilibrium equations are Included. NASTRAN actually uses 
in the stresses at the end of the load step, requiring Iteration to obtain 
this solution. The solution obtained Is a correct nonlinear solution pro- 
vided only that the nonlinear strain contributions to the final stress state 
are negligible. The presence of a dominant Initial stress state Is not needed, 
though It would increase accuracy with this sort of approach. A variant of 
the NASTRAN procedure Is a stepwise calculation in which the stresses at the 
start of each Increment are used In K^. This would eliminate or reduce the 
need for iteration, at the expense of accumulating error* 

Another simple approach which has found widespread use In nonlinear finite 
element work also works without reference to any nonlinear strain descr lo- 
tion. In the derivation of the equations of a particular finite element, 
frequently the geometry of the element Is described In a fairly general 
way from which simpler geometries are available by specialization. For 
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I example, a general quadr i later la 1 element might be considered, and Its 

stiffness matrices and the necessary transformations derived, such that the 
element can be considered In an arbitrary shape and orientation with respect 
to Is coordinate system. As another example, a curved beam element arbi" 
trarlly oriented In a reference coordinate system could be derived. In step- 
wise calculations using such generic elements the nodal and shape descrip- 
tions of the elements can be updated after each step, thereby dealing with a 
new set of elements for each pew solution step. By suitably rotating the 
element local coordinate systems, the strain and stress Increments can be 
referred to the same material fibers for each step and accumulated by sum- 
ming indrem’ents. This again' Is a case in which the nonl inear terms in the 
strain equations are retained only to determine the K„ matrix. Effectively, 
however j nonlinear strain behavior is obtained through the generic element 
updating, though of course, with accumulating error due to the stepwise 
linearization. It Is noted that In general total strains and stresses are 
not obtainable with this method except through summing Increments, due to 
the lack of a governing total strain equation. Thus, corrections for accum- 
) ulating error cannot In general be made with this method, and this Is a 

serious fault. It is also noted that the method Is poorly adapted to large 
strain problems. Many nonlinear finite element programs, particularly those 
developed prior to the last five or ten years, have used this approach. It 
has been referred to as an Euler Ian approach and as an updated Lagrangian 
approach, and refinements have been made to eliminate Its Inherent Inaccu- 
racies (Ref. 105). 

If the finite element derivation retains the nonlinear strain equations not 
only for determining K^, but also for determination of all strain, and hence 
stress, data, then It is -said to be based on a total or Incremental Lagrangian 
approach. There are many options for these purely Lagrangian approaches. 

Most approaches linearize the Incremental strain, by dropping the second degree 
terms in the Incremental strain equations, but retain access to the total 
strain as given by the npnl Inear strain displacement equation. The result 
is similar In the individual steps to the generic element updating scheme 
described above, but wit.h Important advantages and disadvantages. An Important 
I advantage is that the availability of the total Lagrangian strain equations, 
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referred to the same coordinate systems as the Incremental ones, makes It 
feasible to compute exact total strains, thereby avoiding the cumulative 
errors of summing linearized Incremental values. This option Is usually 
exercised by making periodic equilibrium corrections at the end of the com- 
puted Increments. To accomplish this, total stral-hs are used to compute 
element stresses and loads, and from these data residual, or error, loads 
are determined to be applied prior to, or during, the next Increment. An 
Important disadvantage of the total or Incremental Lagranglan procedure Is 
that It suffers from the role exchange of the displacement components which 
was described earlier and Illustrated by Figure 1 . It Is noted that the 
generic element updating method avoids this difficulty due to having 
sequentially rotated the element coordinate system during the stepwise 
solution calculations. 

Procedures such as the purely Lagranglan one. In which fully nonlinear strain 
equations are available, permit the development of. 'finite element formula- 
tions In which the individual Increments are not linearized. In this case 
the Incremental force versus deflection equations are of quadratic (or higher) 
degree, and special solution procedures must be used to solve problems. Vos 
^Ref . 108) has developed such a procedure for a second degree nonlinear step 
behavior. 

One further approach to nonlinear finite element formulation Is discussed here 
This approach combines the advantages of the generic element updating and 
the purely Lagranglan- incrementation. It was developed to avoid the problems 
of displacement component role exchange and stepwise accumulating error, and 
It permits the use of simplified nonl Inear strain equations adapted to beam, 
plate, and shell problems. This simplification, which need not necessarily 
be used, I s the omission of all nonlinear strain terms except those associated 
with lateraV' (bending type) displacements. The procedure Involves the updat- 
ing of the element local coordinate system and transformation of element 
accumulated displacements, but not updating the generic form of the element. 
Thus, the accumulating deformations of the element are followed In a Lagrang- 
lan sense but are referred to a new element local system to start each step. 
Special coordinate transformation equations are required to transform the de- 
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formations in this way, and through these transformations the calculation of 
exact, total Lagrangian strain is accomplished. This allows equilibrium 
checks and corrections to avoid errors due to accumulating linearized in- 
crements. The procedure is described in Volume II of this report. 

All of the approaches briefly discussed above, and others as well, make use 
of the nonlinear strain equations in some way, but differ from one another 
in precisely how the nonlinear terms are incorporated into the calculations. 
Thus the form of the nonlinear strain-displacement equations used in the 
derivation of an element, and also the stiffness matrices of the element, 
are not definitive concerning the true nonlinear character of the element. 

To ascertain this, it is necessary to determine the details of the stepwise 
or iterative calculation procedure used with the element, and to evaluate 
carefully how the approximat ions In both the element der i vati on and the solu- 
tion procedure affect its behavior. 
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6.0 ELEMENTS 



HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-1 

ELEMENT TYPE: Siringer or Truss Memb\^r 


ASSUMED DISPLACEMENT SHAPE: 

' U.= (Xo 
V - ba i' b,K 
Co-^C,X 


Z,0Lr 



DESCRIPTION: 

This element is an axial load carrying member with 2 nodes and 
3 DOF (u, V; w) per node. A stretching deformation state is 
accounted for. Only geometric non linearities are considered. 


REFERENCE: [1], [2], J21] , .[64], 166], [69] . 


VARIATION OF THIS ELEMENT: [92], 193]. 
ADVANTAGES OR DISADVANTAGES: 
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HANDBOOK (CONTINUEDf 


ELEMENT ID: B-T 



DISCUSSION: 


This element is the simplest of the finite elements. It has been derived 
a number of times with fev/ variations, A linear stress-strain relation is 
assumed and only quadratic terms in displacement gradients 'are retained 
in the strain energy. A Taylor series derivation is employed in Ref. [1]. 
The stiffness matrix obtained is termed the initial stress stiffness matrix 
and is reproduced below. A Lagrangian derivation is employed along with 
an incremental solution procedure. 



M 

L 


I -/ 
-/ / 


UL\ 





1 -I o o 

fv.: 

- . - 
u 

f« -7 

- 1 ’ I o o 

0 o 1 -) 

/W,v 



0 O -1 1 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-2 

ELEMENT TYPE: Straight Beam 


a v“ 

► 

Ol 


b 


ASSUMED DISPLACEMENT SHAPE: 





DESCRIPTION: 

This element is a straight beam of constant cross -section, with 2 nodes 
and 3 DOF (u, v, & ) per node. Stretching and bending in the x-z 
plane are accounted for (beading In the x-y plane is an elementary ex 
tension). Only geometric noniinearities are considered. 


REFERENCE: [21, [3], [4], [5], [15], PSl, .164], [66], 168], [88], 

[92], [96] . 

VARIATION OF THIS ELEMENT: [111, [U] 


ADVANTAGES OR DISADVANTAGES: 

Ref [6],".... unsatisfactory for predicting post-buckling behavior 
although adequate for predicting buckling loads." 
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HANDBOOK (CONTINUED) 


ELEMENT ID: B-2 

STRAIN DISPLACEMENT EQUATIONS: 

f- — 


DISCUSSION: 

This element employs a linear stress-strain relation with only quadratic 
products of displacement gradients retained in the strain energy, A 
Lagrangian derivation is employed along with an incremental solution 
procedure. The Initial stress stiffness matrix is reproduced below. 



30L 


36 

3L 

‘36 

3L 


Hi} 



In Reference [3] this matrix is used in an eigenvalue buckling solution 
for beams of variable cross-section. In Reference [4] it is employed vdth 
the midpoint tangent incremental approach, and in Reference [5] if is 
used with nonhnearly elastic materials. In Reference [18] this matrix h 
called the stability coefficient matrix, Reference [64] the nonlinear 
coupling between bending and torsion is included in the geometric 
stiffness matrix derivation. In addition, elastic, ideally plastic, 
linear strain hardening, or nonlinear strain hardening stress-strain 
laws may be used with the element. Viscoelastic effects are included 
in Ref . [68] . . 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-2a 

element TYPE; Straight Beam 


REFERENCE: [11], [37] 





7 


, u 


X, u. 


VARIATIONS FROM BASIC ELEMENT: 

All combinations of displacement gradients are retalnted in the strain 
energy yielding the K, N1 and N2 matrices. Various solution 
procedures are derived and examined in Reference [111. Initial imper- 
fections and an energy perturbation solution approach are used In 
Reference [37], Only geometric nonlinearities are considered. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 

ELEMENT ID: B-2b 

ELEMENT TYPE; Straight Beam 


REFERENCE: [14] 

VARIATIONS FROM BASIC ELEMENT; 

.This is a derivation using convected coordinates in the manner of Argyris. 
It is derived for large displacement, smqli strain problems with material 
nonlinearities. It uses a direct nodal force computational scheme. 

Uses an unusual approach with convected coordinates which separates 
rigid body and deformational displacements through transformation. 
Lagrangian, Eulerlan, and convected coordinates are all used. It is 
a total defornration rather than incremental deformation derivation, and 
it is fully nonlinear while still having linear strain-displacement equations 
Uses direct, numerical integration node force evaluation. Does not make 
use of "initial stresses"; they do not appear in this derivation 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


element ID: B“2c 

ELEMENT TYPE: Straight Beam 



REFERENCE: [63] 

VARIATIONS FROM BASIC ELEMENT: 

This derivation employs the concept of natural loads and natural 
stiffness to formulate the large displacement problem. It is 
derived for large displacement, small strain problems using con- 
vected coordinates. An incremental solution procedure is developed 
The method is applicable to post-buckling analyses where an itera- 
tive procedure is used to find the initial post-buckled position. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-2d 

ELEMENT TYPE: Straight Beam 


REFERENCE: Il20] 



VARIATIONS FROM BASIC ELEMENT: 

■L dv 2 1 du 2 

Author retains both *2 ^ dx ^ *2 ^ dx ^ strain terms, though the latter 

proves unimportant in his examples. Principal variation from the basic element 

is that the loads p(x) and P ,P are permitted to change direction as functions 

X y 

of the deformation, producing ’’follower load" effects. With these effects re- 
presented by a matrix L, authors formulation has the form 

AQ = [[K] + [G] + {LJJ Aq 

where G is the geometric stiffness matrix and K is the stiffness matrix for 
linear behavior. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-<3 

ELEMENT TYPE; Straight Beam 


ASSUMED DISPLACEMENT SHAPE: 


u,v 



V - h,X 


DESCRIPTION: 

This is the "stability" beam element where the quintic u displacement 
is consistent with the cubic v displacement from the equation of axial 
equilibrium. It has two nodes and 3 DOF (u, v, v ) per node. Stretch- 
ing and bending deformations are accounted for. Only geometric non- 
linearties are considered. The excess generalized DOF are reduced on 
the elemental level. 

REFERENCE; [6] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
Improved postbuckling behavior. 
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DISCUSSION: 

This element employs a linear stress-strain relation with all displace- 
ment gradient combinations (up to quartic) retained in the strain energy. 
A Lagrangtan derivation is employed. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-4 

ELEMENT TYPE: Curved Beam 



ASSUMED DISPLACEMENT SHAPE: 





-S \Vn(j 


0 

0 


tpS (4^^) C) 




At 


i- U 




O o 
yi ~ip3 

zy 3Y 


C<^:\ 

I 

aj 

0.5" 




DESCRIPTION; 

This is a curved beam with 2 nodes and 3 DOF (u, w, ^ ) per node. 

The rigid body modes are added explicitiy to the strain -^inducing modes. 
Stretching and bending deformation states are accounted for. Geometric 
and material nonlinearities are considered. 


REFERENCE: [7], 112] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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HANDBOOK (CONTINUED) 


ELEMENT ID: B-4 

STRAIN DISPLACEMENT EQUATIONS: 



DISCUSSION: 

The element is developed for an elastic-plastic, strain -hardening, and 
strain -rate material behavior through the use of the mechanical sublayer 
model. The element is developed for a transient response analysis and 
the form of the dynamic equations are not standard. Matrices are 
developed in terms of the products of stresses with linear and nonlinear 
strains to obtain well behaved numerical integration. A Lagrangian 
derivation is employed and all nonlinear products are retained in the 
Virtual y/ork expression. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 

ELEMENT ID: B“5 

ELEMENT TYPE: Circular Arch 



ASSUMED DISPLACEMENT SHAPE: 

fv 1 eos ("f-) ^s\n('k) 1 ^ 

\oj) _<i>> ^ (f') 0 R 



DESCRIPTION: 



per node. A constant circumferential strain and a linearly ' 

varying bending strain is assumed. The deformation , states arise from 
the simultaneous solution of these equations. Only geometric non- 
linearities are considered. 

REFERENCE: [8], [9] 


VARIATION OF THIS ELEMENT: 

ADVANTAGES OR DISADVANTAGES; 
improved convergence. 

Probably not suitable for large element deflections or for 
post -buck ling/ due to constant circumferential strain 
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HANDBOOK (CONTINUED) 




ELEMENT ID: B-5 

STRAIN DISPLACEMENT EQUATIONS: 

^ cJi/ uj* f / \ 

'^isjy "4 ( ip- ' K- y 


DISCUSSION: 

In Reference [8] the linear element of Reference [9] is extended to 
geometrically nonlinear problems. Rigid body states are explicitly formed 
because of the integration of the assumed strain states. The rational 
for the choice of strain states is apparently empirical. All combinations 
of displacement derivatives are retained in the strain energy. The 
Lagrangian strains are employed. The element is used in the linearized 
incremental method based on mid-increment stiffness and Newton -Raphe on 
iteration. 



HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID; B-6 

ELEMENT TYPE; Shallow Curved Beam 


/ 



ASSUMED DISPLACEMENT SHAPE; 

W= [ )j X, 0,0] 

~ Not specified, needs to be dealt with in matrix generation, 
by integrati6ns. 


DESCRIPTION; 

This is a shallow curved beam with 2 nodes and 4 DOF per node 
(u, du/dx, w, dw/dx). An initial imperfection, w , is used to 
specify the initial curved geometry. Stretching anc? bending deforma- 
tion states are accounted for. Only geometric nonlinearities are 
considered. 

REFERENCE; [ 10 ] 


VARIATION OF THIS ELEMENT; H2] 


ADVANTAGES OR DISADVANTAGES; 
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HANDBOOK (CONTINUED) 


ELEMENT ID: B-6 

STRAIN DISPLACEMENT EQUATIONS: 



(not explicitly presented) 
(potential energy method used) 


DISCUSSION: 

This element employs a linear stress-strain relation with only the 
quadratic combinations of displacement gradients retained in the 
strain energy, thus yielding the classical k' matrix. It should be 
noted that retaining the displacement gradient du/dx as a nodal 
DOF does not conform with the usual element assemblage procedure. 
A Lagrangian derivation is employed. The mid-increment solution 
procedure is proposed. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-6a 


ELEMENT TYPE: 


BEAM (For in-plane large displacement analysis) 

^/ern-<L>'\i 0. c.r 



where ? = z/h 

0 ; left end 

1 ; right end 


DESCRIPTION: 2 nodes, 4 degrees of freedom per node - 

(Stretching and bending deformations are modelled) . 

Displacement components are referred to the initial 
configuration. 


31 



HANDBOOK (CONTINUED) 


ELE3NT ID; 

ADVANTAGES OR DISADVANTAGES: The cubic polynomial representation for 

both displacement components, allows for the use of this element with 
arbitrarily large displacements and rotations. (The use of a linear 
representation for the 2 - component would lead, on the contrary, to a non-* 
obj active curvature measure ruling out the use of such an element for 
cases of large rotations). 

STRAIN DISPLACEMENT EQUATIONS: 

e = w' + -^ (w'^ + u'^) 

X = u”(l + w') - w” u* 
where e: axial stretching measure 

X: axial bending (curvature) measure 

DISCUSSION: The element is used in connection x^ith an exact non-linear 

formulation analogous to Budiansky’s non-linear Shell Theory (ASME T. of 
Appl. Mech. , June 1968). The strain-energy expression used is equivalent 
to the 3-D expression under Kirchhoff’s hypothesis. Constitutive equations 
may be prescribed in terms of Kirchhoff’s stress tensor and Green’s strain 
tensor and integrated over the cross-section. 

Resulting equations are polynomials of the 3rd degree in 
the nodal displacements and their first z-derivatives. Preferred solution 
procedure: Newton-Raphson. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-7 

ELEMENT TYPE: Curved Beam Element 


ASSUMED DISPLACEMENT SHAPE: 




Cos (f) 





o 

V o o 


DESCRIPTION: 

•V y 

This is a curved beam with 2 nodes and 4 DOF ( V, W, r- ’‘%-9 -X, »nd 

) per node. The rigid body modes are added explicitly 
to the strain inducing modes. Stretching and bending deformation states 
are accounted for, i.e., a Timoshenko type element. Geometric and 
material nonlinearities are considered. 

REFERENCE: [12] 

VARIATION OF THIS ELEMENT: HO] 

ADVANTAGES OR DISADVANTAGES: 

Improved strain prediction. 



HANDBOOK (CONTINUED) 




ELEMENT ID: B-7 

STRAIN DISPLACEMENT EOUATIONS: 





V_\^ 
R f 


DISCUSSION: 

The element Is developed for an elastic plastic, strain hardening, and 
strain rate dependent material behavior through the use of the mechanical 
sublayer model. The element is developed for a transient response analysis. 
All nonlinear products are retained in the virtual work expression. The 
author states the fourth DOF, c^~ , at each node is retained because of a 
lack of a rational condensation process for dynamic problems. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B~8 

ELEMENT TYPE: Curved Beam Element 


ASSUMED DISPLACEMENT SHAPE: 






~2 do6C^ ^ 


t 

0 

o 

o 

J. 


<j> 


o 

o 

Y 

o 


TJ 

0 

0 

O 

o 

o 

0 

1 

? 


DESCRIPTION; 

This is a curved beam with 3 nodes and 3 DOF (v, w, S) at the two 
end nodes and 2 DOF (v, w) at the midpoint node. The rigid body 
modes are added explicitly to the strain inducing modes. Stretching, 
bending, and shearing deformation states are accounted for. Geometric 
and material nonlinearities are considered. 

REFERENCE; [12] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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HANDBOOK (CONTINUED) 


ELEMENT ID: B-8 

STRAIN DISPLACEMENT EQUATIONS: 

C - 4.^ ^ 4. ' ^ 

ly - y r/4^ U ^ 


DISCUSSION: 

The element is developed for an elastic-plosHc, strain hardening, and 
strain rate dependent material behavior through the use of the mechanical 
sublayer model. The element is developed for a transient response analysis^ 
All nonlinear products are retained in the virtual work expression. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-9 

ELEMENT TYPE: Curved Beam Elemenf 


ASSUMED DISPLACEMENT SHAPE: 




"H. dCC' ^ f I’l ^ O O 

O O 

o o i ^ 



DESCRIPTION: 

This is a curved beam with 2 nodes and 3 DOF (v, w, 0) per node. 
The rigid body modes are added explicitly to the strain inducing 
modes. Stretching, bending, and shearing deformation states are 
accounted for. Geometric and material nonlinearities are considered. 


REFERENCE: [12] 


VARIATION OF THIS ELEMENT: [i6] 


ADVANTAGES OR DISADVANTAGES: 

Requires small element size to accurately obtain load output. 
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HANDBOOK (CONTINUBP) 


ELEMENT ID: B-9 


STRAIN DISPLACEMENT EQUATIONS: 


£ <7 J) •£.<?)* , r(7,:f)-K(7) 

K-lf , 


DISCUSSION: 

The element is developed for an elastic-plastic, strain-hardening, 
and strain rate dependent mdterial behavior through the use of the 
mechanical sublayer model. The element is developed for a transient 
response analysis. All nonlinear prodvcts are retained in the virtual 
v/ork expression. 


38 



HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-9a 

ELEMENT TYPE: Straight Beam 


REFERENCE: [16] 



VARIATIONS FROM BASIC ELEMENT: B“9 

This is also a mechanical sublayer model, but the element Is generated 
from an assembly of flanges and shear webs and not from energy consid 
erations. The element is developed for an elastic - perfectly plastic 
material and a transient response analysis^ Geometric nonlinearities 
are accounted for by updating coordinates. 
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REPBODUCIBILITY OF' 
OIKMAL PAGE IS FO'.:-: 


HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-10 

ELEMENT TYPE: Stringer or Truss Member 


ASSUMED DISPLACEMENT SHAPE; 



u(x)' 

v(x) 

. = (2^^ - 3c + 1) < 

[" 1 
Y. 

1 

4- (2c^ " zy 

[ n . ' 
3 

+ 4(5- 5^) , 


. w(x) 


w . . 
X 

i 

Iw. , 
3 


[wkl 


Initial strain distribution e 


= constant 


DESCRIPTION: 

This element is a linear strain, axial load carrying member with 3 
nodes and 3 DOF (u,v,w) per node. A stretching deformation state 
is accounted for. Both geometric and material nonlinearities are 
considered. 


REFERENCE: [64] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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HANDBOOK (CONTINUED) 


ELEMENT ID: B-10 

STRAIN DISPLACEMENT EQUATIONS: 

= ^ + i 4 . i 

3x 2 i QxJ 2 V 3x/ 


DISCUSSION; 

This element employs a linear elastic stress-strain relation with 
only quadratic products of displacement gradients retained in the strain 
energy. A Lagrangian derivation is employed to yield the geometric 
stiffness matrix [k^] . The ^tiffness and geometric stiffness matrices 
are reproduced below (where F is the average axial force) . Apparently 
the updated' or convected coordinate approach to the solution of the 
geometrically nonlinear problem is used. Plastic strains are cal- 
culated using a uniaxial bilinear stress-strain relation or a nonlinear 
Ramberg-Osgood relation. Perfect plasticity is also accommodated. 


[k“J{6} = 



[k’]{5> = 



1-8 
7-8 
8 16 
0 0 

0 0 

0 0 - 


0 0 0 

0 0 0 

0 0 0 

7 1 - -8 

1 7 - 8 

8-8 16 



A1 
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HANDBOOK OF NONLINEAR 
FINITE ELEMENTS 


ELE::EHT ID; B' 11 

ELEMENT TYPE; Beam (Linear elastic) 

ASSUMED DISPLACEMENT SHAPE: 

1. Cubic polynomial for u displacements. 

2. Cubic polynomial for v displacements. 

3. Cubic polynomial for torsional displacements. 

4. Linear polynomial for axial displacements. 

5. Cubic polynomial for flange displacements. 

DESCRIPTION: 

I 

(a) For member buckling the element has 14 degrees of freedom 
associated with the nodal variables shown in Fig. 1. 



< y > i < Vi , iv.il , Vj , > 

< > i < fp, , , ‘f^j , ^ 

< W > : < Wf , V<' J > 


Figure 1 - Beam Element and Nodal Displacements 
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(b) For local buckling the element has 8 degrees of freedom, 
four associated with the top flange (as illustrated in 
Fig. 2), and four associated with the bottom flange. 



< I > = < < > > 


Fi gure 2 - Plate Flexure Shape Functions and Nodal Displacements 


The element is formulated with respect to the initial coordinate system 
(Lag ran gi an formulation). 


REFERENCE: (7f)j 


VARIATIONS OF THIS ELEMENT: 

Considering only the (a) displacement patterns the element can 
be used for a linear eigenvalue member stability analysis which couples 
axial, torsional and biaxial bending effects. Introducing the (b) 
displacement pattern allows local flange buckling to be coupled to the 
behavior in (a). When (b) displacement patterns only are used, an uncoupled 
flange buckling analysis may be performed. 


3 . 




ADVANTAGES OR DISADVANTAGES: 

Beam buckling problems or load deformation problems Involving 
three dimensional response may be modelled with a minimum number of 
degrees of freedom. 

STRAIN DISPLACEMENT EQUATIONS: 

(a) For member displacements, 
u = Uq — (y - e^) tj) 

V - v^ + (k - e^) <|) 

w = - y v; - X u • + 4) 

. _ BW 

BZ 

in which, and v^ are displacements of the shear center; 

VI ^ is the displacement of the centroid (origin of x - y 
coordinates); and 4) is the twist. 

(b) For local plate displacements, 

u* = - zwt 

’a 

v^ = - zwto 

p 

in which w'*' is the traverse local plate displacement, 
a and 3 are local coordinate axes, and z is the normal 
coordinate measured from the mid-surface of the plate. 


DISCUSSION: 

Stress -strain response is assumed linear elastic. Derivation 
is by the principle of virtual displacements with the majority of terms 
in the "geometric stiffness” matrix being retained. Solution may be 
either a linear eigenvalue solution or a nonlinear load-deformation 
solution. Lagrangian coordinates are used throughout. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT IE: B*>12 

ELEMENT TYPE: Inelastic beam (Shanley tangent modulus) 

ASSUMED DISPLACEMENT SHAPE: 

1. Cubic polynomial for u displacement. 

2. Cubic polynomial for v displacement. 

3. Cubic polynomial for torsional displacement. 

4. Quadratric polynomial for axial displacement. 


DESCRIPTION: 

The element has two nodes at each end, one for u, v and cj) 
nodal degrees of freedom, and one for a w degree of freedom (see Figure 1) 
In addition there is a midpoint node for an additional w degree of freedom 
which permits a quadratic polynomial to be used''for axial displacements. 
This degree of polynomial is the minimum permissable to obtain a balance 
of nodal forces at interface nodes when inelastic deformations occur. 

The element tangent stiffness has been formulated by employing 
transformed section properties. A tri linear stress strain curve, as 
shown on Fig. 2, allows subdivision of each component plate segment into 
at most five regions, as shown in Fig. 3, which may be assembled into a 
section of equivalent stiffness as shown in Fig. 4. 



it) - < l^x > 

(w) = < w§ > 

Figure 1 - Finite Element Degrees- 
and Nodal Vectors 



Figure 2 - Tril inear Stress-Strain 
Diagram 


-Freedom 
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2. 



Figure 3 - Transformed Section of Figure 4 - Arbitrary Transformed 
Plate Segment Section 

REFERENCE: I71], 


ADVANTAGES OR DISADVANTAGES: 

■^fie formulation refers element properties to a fixed set of 
axes so that incremental element matrices do not require a transformation 

prior to assembly. 


strain DISPLACEMENT EQUATIONS: 
u = Uq - (y - e^) (!) 

V = + (x - e ) ({) 

0 X 



in which Uq and are displacements of an arbitrary 'shear centerVaxis 
w is the displac^ent of an arbitrary ‘centroid axis' (origin of x - y 
coordinates); and (j) is the twist of the cross-section. 
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3 . 


DISCUSSION: 

Derivation of equilibrium equations is by the principle of 
virtual work and retains terms significant in the geometric stiffness 
matrix. An iterative Newton Raphson procedure is used for solution but 
equilibrium is checked against the total applied load, thus preventing 
drift from the true solution. The program can be used for an eigenvalue 
solution to give inelastic buckling loads or for a nonlinear load 
deflection analysis to give inelastic instability loads. Lagrangian 
coordinates are used throughout. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-I3 

ELEMENT TYPE: FLA 2 (normal stress flange In 3“space) 



ASSUMED DISPLACEMENT SHAPE: 1 i nea r 


DESCRIPTION: 

Number of nodes : 2 

Degrees of freedom : u, v, w at each node, all together 6 

Deformation type : stretching 

Isotropic 


REFERENCE: 1^2] , l93l 


VARIATION OF THIS ELEMENT; IlJ, l2] , l2l], I6A] , 166], 169], 


ADVANTAGES OR DISADVANTAGES: 

DISCUSSION: Based on the concept of natural stress invariants for large 

rotations, 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-1^ 

ELEMENTTYPE: Elastic rod 


ASSUMED DISPLACEMENT SHAPE: Plece~wlse linear; contains built 

in discontinuities for shock wave calculations, • 


DESCRIPTION: 

No. of nodes: 2 

Degrees of freedom; 4 
Coordinate system: Cartesian 

Isotropic materials 


REFERENCE: 

iiio], iiii],'ni2], 1113], 1114 ], 



VARIATION OF THIS ELEMENT: 

ADVANTAGES OR DISADVANTAGES: 

Allows for sharp definition of shock fronts without an undesirable degree 
of artificial dissipation. 
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strain Displacement Equations : 




j 


+ u . . + u . u . ) 

3 , 1 . m,i 


Discussion t 

1. Stress-Strain equations: ■) + 

2. Terms retained in the strain energy: All 

3. Derivations from other than strain energy considerations: 

Gaierkin Methods, Conservation of Energy. 

' . ^ 

4. Unique derivation techniques: See (c) above. 

5. Lagrangian, Eulerian, convected: Lagrangian (convected on 

material) 

6. Numerical integration schemes: 

For static problems: Incremental Loading (Euler's Method) 

plus Newton-Raphson corrections. 

7. Required or preferred solution procedure: 

For dynamic problems; Incremental Loading. Divided central 
difference methods plus a Lax-Wendroff method. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-15 

ELEMENT TYPE: bECOS (BECOSX) 

(solid beam of constant cross section in 3 space) 



DESCRIPTION: 

Number of nodes 
Degrees of freedom 

Deformation type 

Isotropic 


located 
V , w. 


2 , 

U, V, W, (j) , 
tpgether 1 T 


on centroid of end cross-sections 
^ , c{) at each node, all 


stretching, bending without shear effects, 
St. Venaht Torsion 


REFERENCE: [92], [96]. 


VARIATION OF THIS ELEMENT: [2], [3j , [4], [5] , ‘ [15] , [18j, [2lJ, [64], [66], 

[68], [88], The concept was also applied to beams connected to eccentric 
nodal points (BECOSX). 

ADVANTAGES OR DISADVANTAGES: 

I 

Strain-Displacement Equation and Discussion : The original derivation of this 

element is based on the natural stress invariants for large rotations. One ' 
special version started also with nonlinear strain-displacement relations 
(large rotations;. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: B-16 

ELEMENT TYPE: CIRCA 

Ccircular beam in space) 


ASSUMED DISPLACEMENT SHAPE: Third order 

Figure ; 


2 

V, w, (j)^, (J)^, at any nodal 

point, all together 12 
stretching, bending torsion (St. Venant) . 
. \ 

REFERENCE: |g6] 

VARIATION OF THIS ELEMENT: 

ADVANTAGES OR DISADVANTAGES: 

D I SCUSS I ON : Based on the concept of natural stress invariants for large 

rotations. 


DESCRIPTION: 

Number of nodes : 
Degrees of freedom: 

Deformation type 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: 
ELEMENT TYPE: 



‘ ASSUMED DISPLACEMENT SHAPE: 

Ui- a,x ^ 

v= *h,K. , ■ 

6r = C„*C,y+Qy • 

DESCRIPTION: 

This element Is a plane stress triangular plate with 3 nodes and 3 DOF 
(u, V, w) per node. Only the stretching deformation state is accounted 
for. Only geometric nonlinearlties* are considered. 

REFERENCE: HI, (2], I21J,I60J 

1.107], [113], in'll, Ills], 1117], [ 118 ], t9l|]. 

VARIATION OF THIS ELEMENT: 

ADVANTAGES OR DISADVANTAGES: 


Fully compatible element 





HANDBOOK (CONTINUED) 

ELEMENT ID: P-1 


STRAIN DISPLACEMENT EQUATIONS: 



3<xJ 


DISCUSSION: 

This is another of the simplex elements with an assumed linear dis- 
placement state which results in a constant stress, constant strain element. 
A linear stress -strain relation is assumed and only quadratic terms are 
retained in the strain energy. A Lagrangian derivation is employed along 
with an incremental solution procedure. Node force increments are 
determined from stiffness matrices. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-la 

ELEMENTTYPE: Membrane and Bending 

Triangular Plate 


REFERENCE: [20], [ 21 ] 

VARIATIONS FROM BASIC ELEMENT: P“1 

This element includes bending and shearing deformations in addition to 
the stretching deformation states described in P-1, The additional 
assumed displacement states are: 

6 ^- e, X + 

where 

u.= uU^ij) 

V = G.y) 

Only geometric nonlinearities are accounted for. The derivation is 
generalized to Include orthotropic materials. 

Degrees of freedom per node are u, v, w, 0 , 0 , 

X y 




HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-lb 


ELEMENT TYPE: Membrane Triangular Plate 


REFERENCE: [ i 9 ] 



VARIATIONS FROM BASIC ELEMENT: P-l 

This derivation considers material nonlinearities in conjunction 
with the basic isotropic membrane plate. Strains are assumed to 
be of the form 

= T 

Ctj^ ' oil 

where the subscripts e and pt denote elastic and inelastic 
respectively. The inelastic strains include both creep and time- 
independent plastic strains. The inelastic strains are converted 
to a column matrix of loads. An incremental solution procedure 
is employed. 
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handbook of nonlinear finite elements 


ELEMENT ID: P-lc 

ELEMENT TYPE: Membrane Triangular Plate 


REFERENCE: [U] 

VARIATIONS FROM BASIC ELEMENT: P-1 


A 



This derivation considers large rotations in the plane of the element. 
Convected coordinates are used which separate rigid body and de- 
formation displacements through transformations. Lagrangian, Eulerian, 
and convected coordinates are all used. It is a total deforrrKition 
rather than an incremental deformation derivation, and it is fully non- 
linear, yet having linear strain -displacement equations. It uses direct 
numerical integration node force evaluation„ It is derived for large 
displacement, small strain problems with material nonlinearities. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-2 

ELEMENT TYPE; Rectangular Plate 


z>os 


ASSUMED DISPLACEMENT SHAPE: 

^ f— CiJ 



C2.,Z) 


•j^ X, U- 


Oj(x,ij)= ^ H,jCy)4Ty 

A/oiCy.1 ~ ^2.)(^~3a5^^ta^)/a^ H,!' U)~ /a^ 

= ^K^-ax^a" 

Similarly for y direction replace x by y, and a by b. Similar 
expressions for u and v. 

DESCRIPTION: 

This element is an isotropic rectangular plate with 4 nodes and 12 DOF, 

(u, V, w, u^, v^, w^, u^, V , Wy, u , Vj^y, w y) per node. Stretch- 
ing and bending deformation sfates are accounted tor. Displacements 
are approximated by products of one dimensional cubic Hermite interpola- 
tion polynomials. Only geometric nonlinearities are considered. 


REFERENCE: 


16] 


VARIATION OF THIS ELEMENT; See also P-12, P-2a, S-5 


ADVANTAGES OR DISADVANTAGES: 

The bicubic interpolation functions for u,v, and w admit to six 
linearly independent displacement states of very little strain energy. 
Consequently, this element offers a better resolution of the rigid body 
mode problem. The addition of the w^y terms to the transverse deflec- 
tion corrects the deficiency of the Papperfuss element [28] in allowing q 
constant twist state. The selected nodal DOF may not be compatible with 
other elements in a general purpose structural analysis computer program. 
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HANDBOOK (CONTINUED) 


ELEMENT ID: P-2 

STRAIN DISPLACEMENT EQUATIONS: 




' c)uf \ 2- 

> 

^ d'^ur\ 

c)^ a ( 

TT / 

2 ( 

) 

S,- 


• dw 


( d^cj\ 


c).jJ - 

2 ( 

^ Tp- J 

K 

^ d u. dv 

dcU" 


c)j ^ da 

■ ^ dx 


- 


DISCUSSION: 

This element employs a linear stress-strain relation with all combina- 
tions of displacement gradients retained in the strain energy. A 
Lagrangian derivation is employed with a numerical solution obtained 
by direct minimization of the total potential energy, 

u, V, w, w , w are all continuous at nodes. Since w is cubic 

, X y X 

in y on an X = constant edge, the specification of w cannot control 

over the entire edge, so that the slopes are not everywhere continu- 
ous (e.g., middle of a side). 

Specification of continuity u,u,v,v,u ,v , and w really 

X y X y xy xy xy 

constitute a stress or strain constraint, which is undesirable, and the element 
could be referred to as "over-constrained", for this reason. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-2a 

element TYPE: Rectangular Plate 


REFERENCE: [38], [39], [40] 

VARIATIONS FROM BASIC ELEMENT: P~2 

This element is a rectangular plate in which, for Ref, [38], complete 
generality is maintained in plate material properties so that anistropy 
or isotropy is considered, as well as plastic yielding prior to buckling. 

It has 4 nodes and 4 DOF (w, w , w , w ) per node. Only bend- 
ing deformation states are considered since ^Ke analysis is limited to 
buckling solutions. The transverse displacement function is that of 
Bogner, Fox, Scbmit, Ref. [25]. The conventional initial stress stiffness 
matrix is developed where only quadratic terms are retained in the strain 
energy. In addition, the standard beam element B-2 is used in conjunc- 
tion with the plate to perform stiffened plate buckling analyses. 

In Ref, [39], the initial stress’ stiffness matrix for the BFS element. 

Ref. [25] is suggested. 

In Ref. [40] the initial stress stiffness matrix for the BFS element is 
developed and applied to the plastic buckling analysis of plates. 



HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: 

ELEMENT TYPE: Rectangular Plate 


ASSUMED DISPLACEMENT SHAPE: 







^)'*’ coi^c^c,: Xp~ o. ^ x^~ 0. 


DESCRIPTION: 

This element is a rectangular plate with 4 nodes and 3 ^pP (^r '^x z- '^y) 
per node. The linear stiffness matrix was first presented by Melosh. 

The matrices are derived for a buckling analysis consequently only bending 

deformation states are required. 


REFERENCE: [17], 136] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 


61 


HANDBOOK (CONTINUED) 


ELEMENT ID: P-3 

STRAIN DISPLACEMENT EQUATIONS: 
(Not explicitly given) 



DISCUSSION: 

This element employs a linear stress-strain relation with only 
quadratic combinations of displacement gradients retained in 
the potential energy „ A Lagrangian derivation is employed with 
an eiganvalue analysis performed to find the buckling loads. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P"4 


^,V 



The author speaks of an assumed linear stress distribution for O' x cind 
f and a constant shear stress distribution. A formal derivation was 

u y 

not presented „ 


DESCRIPTION: 

This element is an isotropic rectangular plate with 4 nodes and 2 DOF 
(u, v) per node. It accounts for material nonlinearities but does not 
consider geometric nonlinearities. Only membrane deformation states are 
required. 


REFERENCE: [19] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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HANDBOOK (CONTINUED) 


ELEMENT ID: 


STRAIN DISPLACEMENT EQUATIONS: Given in general as: 


6. H = 6 ze + c<T ^ 6 jE.pT- 


c^tKy 

6 ^ = 


£ 




is the total inelastic strain component 


DISCUSSION: 


This derivation considers material noniinearities for a membrane plate. 
Plastic strains are converted to a column matrix of plastic loads. An 
incremental solution procedure is employed. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-5 

ELEMENT TYPE: Strip Plate Element 



ASSUMED DISPLACEMENT SHAPE: 


V - i 

^ b 



DESCRIPTION: 

This element is a rectangular strip with two line nodes and two DOF 
(translation, rotation) per node. A sinusoidal distribution for w is 
assumed along the length. Only bending deformation states are accounted 
for since the element is used in an eiganvalue buckling analysis. 


REFERENCE: [22], [23] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 



HANDBOOK (CONTINUED) 


ELEMENT ID: P-5 


STRAIN DISPLACEMENT EQUATIONS: 


•7 

c)^ 


Cf.y 


^ -- 2 £ 


> --J 


DISCUSSION: 


This element employs a linear stress -strain relation. Only quadratic 
products of the displacement gradients are retained in the strain energy 
A Lagrangian derivation is employed with the solution for buckling 
loads obtained by an eiganvalue analysis. 
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ELEMENT ID: P-6 


ELEMENT TYPE: 



ASSUMED DISPLACEMENT SHAPE: 

^ (i>) ip) 


(o) 


U= Z ^o\C<) HoJ (uJUi,; ; /4, '"aCX-aj " similarly for v 

'j J 

p -J-. r 


- L-- 
I-I jr/ ^ 


/i\ A ' Lf^ w {'A . , \ I N / 




<:/3 


^0 


CO 


to 




c<> 


Hoi (>0 - /Voi /H;, r cl^x j 

Similarly for y cJireetion replace x by y, and a by b, 


H / Z ~ cj 

DESCRIPTION: 


f 3, . -L iS> t) 

'C)C -C.X / 


This element is an isotropic rectangular plate with initial imperfections. 

It has 4 nodes and 6 DOF (u, v, w, w , w , w ) per node. Stretch- 

. ' X y X y 

ing and bending deformation states are accounted for. Displacements 
are approximated by Lagrangian and Hermitian interpolation polynomials. 
Only geometric non linearities are considered. 


REFERENCE: [24]; see [25] for linear element^ 

VARIATION OF THIS ELEMENT: 

ADVANTAGES OR DISADVANTAGES: 

Assumed displacement state is complete and therefore converges 
monotonicaliy. 

Linearity of u, v probably means not suitable for large deflection - 
post-buckling type of problem. 
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HANDBOOK (CONTINUED) 


ELEMENT ID: P-6 


STRAIN DISPLACEMENT EQUATIONS: 

c)u I \ .y / ^ \ 

Cx “ 2 L J a Wx / ^ ^ Ox"'/' 

^ - 0^ ; r t /c)uji\z, 

L J 

^ ^ ii! ^ ^LuX^Uo) _ 0^ OuJ^ 5 

r)y 0j< Oy 


DISCUSSION; 

This element employs a linear stress-strain relation with all 
combinations of displacement gradients retained in the strain 
energy. Initial imperfections (w ) are included. A Lagrangian 
derivation is employed. The equations for both the direct iterative 
analysis and the step-by-step incremental procedure are presented. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 
ELEMENT ID; P-7 

ELEMENTTYPE: Rectangular plate or shallow shell 

^,Of 

ASSUMED DISPLACEMENT SHAPE: 

=L ^1 y, J /o<J 

v = L-^:X.y.xyJ Ha] 

Cj xy , x", X V, xy^ /y. xy 'j 

DESCRIPTION: 

This element is a rectangular plate or a shallow shell. It has 
4 nodes and 5 DOF (u, v, w, w , w ) per node. Stretching 

\ t • f y 

and bending deformation states are accounted for. Only 
geometric non linearities are considered, 

REFERENCE; [26] 

VARIATION OF THIS ELEMENT: [10], [45] 

ADVANTAGES OR DISADVANTAGES: 

It does not satisfy normal slope compatibility. 
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HANDBOOK (CONTINUED) 


ELEMENT ID: P-7 

STRAIN DISPLACEMENT EQUATIONS: 

^z. " ^ i^y ^ ^ 

(f,= U,,^ -f Ui,i< -2 -Sor.xj, 

DISCUSSION: 

This element employs a linear stress -strain relation with all combinations 
of displacement gradients retained in the strain energy. A Lagrangian 
derivation is employed with the numerical solution obtained by Newton- 
Raphson iteration. The formulation is also specialized to obtain the 
linearized incremental equations and the equations for a linear buckling 
analysis. Initial strains are considered. 
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ELEMENT ID: P-7a 


ELEMENT TYPE: Rectangular Plate or Shallow 



This element is an isotropic rectangular plate with initial imperfec- 
tions or a shallow shell. The displacement shapes are the some as 
element P-7, but only quadratic products of displacement gradients 
are retained in the strain energy. The integrals are evaluated 
numerically using Gaussian quadrature formulae. The solution is 
obtained using mid-increment stiffness matrices. Only geometric 
nonlinearities are considered. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P-7b 

ELEMENT TYPE: Quadrilateral Plate 

I ' 

REFERENCE: [45] 

VARIATIONS FROM BASIC ELEMENT: P"7 



This element uses the membrane stiffness matrix of Ref. [19], (see P-4), 
with the 12 term polynomial of P-7 for the lateral displacement w. A 
linear stress-strain relation is employed with only quadratic products of 
displacement gradients retained in the strain energy. Several assumed 
transverse displacement functions are used in deriving the k' matrix 
including 1) the 12 term polynomial used in developing the kP matrix, 

2) a six term, assumption, and 3) a four term assumption. A Lagrangian 
derivation is employed with the numerical solution obtained by a com- 
bined incremental step and eigenvalue analysis to determine the buckling 
load. 
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ELEMENT ID: P-8 

element TYPE: Rectangular Plate 



ASSUMED DISPLACEMENT SHAPE: 
-V 3 


k~{ -<~l 


DESCRIPTION 


This element is an isotropic rectangular plate with 4 nodes and 3 DOF 
(w, p< f j6 ) per node. Only bending deformation states are accounted 
for. An initial state of constant stress is assumed and the resulting 
initial stress stiffness matrix issuitable for a buckling analysis. 


REFERENCE: [27], [28] 


VARIATION OF THIS ELEMENT; [29 ] 


ADVANTAGES OR DISADVANTAGES: 

This element omits the deformation state corresponding to constant twist 
and although it is fully compatible and will converge monotonical ly, 
it may not converge to the true solution. 


HANDBOOK (CONTINUED) 


ELEMENT ID: P-8 


&TftAHt€>t9Pt:AeeMENT‘eQttATfmS; 
ENERGY EQUATION: 


u- f 


II {^- 






DISCUSSION: 

This Is the Papenfuss element [28], extended to a buckling analysis. 

It employs a linear stress -strain relation with only quadratic combina- 
tions of displacement gradients retained in the work expression. A 
Lagrangian derivation Is employed with an eigenvalue analysis per- 
formed to find the buckling loads. 
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ELEMENT ID: P-8a 

ELEMENT TYPE: Rectangular Plate 


REFERENCE: '[29] 



VARIATIONS FROM BASIC ELEMENT; 


This is the element of Reference [28] extended for use in the 
large deflection of plates, A membrane displacement shape is 
assumed 

V - Ll- -^)v; ^ (i~ 


and used in conjunction with the bending shape of P-8 in the strain 
displacement equations 



57T 


'cly ^ 

vl - ^ f o -T 

()a ^ cVc)y 


The resulting element has 4 nodes and 5 DOF (u, v, w, , /6 ) per 
node. A linear stress-strain relation is used and all combinations of 
displacement gradients are retained in the strain energy. The i)Vird 
and fourth order combinations of displacement gradients are converted 
to equivalent load columns and the resulting equations solved by an 
iterative technique. A Lagrangian derivation is employed. 
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ELEMENT ID: P-9 

ELEMENT TYPE: AnnuloF- Plate * ■ 



ASSUMED DISPLACEMENT SHAPE: 


K = ^ = Vo<V>7ck/A ^ 2('c<-/)/6 

^ Y?c- = i fao<"/nc< -o(^-pi)/A 

" iio{ln cK 0/A , ~ I - h<^{ 2.oi^!in - o< + | ^/l\ ^ Y/‘/ -h<^ ^ 2. Jr) o( )//\ 


DESCRIPTION: 

This element is an isotropic annular plate with 2 ring nodes and 2 DOF 
(translation, rotation) per node. The displacement function is the exact 
solution to the circular plate equation -bending only. An initial state of 
stress is assumed and the resulting initial stress stiffness matrix is suitable 
for a buckling analysis. 

REFERENCE: [30] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 


■ 
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1 ELEMENT ID: P-9 

ST1tAftTB1«PLTXG€MENT-EQU'ATfON^ 

PotenHal Energy: 

V= iJ^2 7T Cot')^ ^d[(0J''+ - 

cC~ f 

- No cC'-'b^ [N I 


DISCUSSION: 

This element employs an isotropic, linear stress-strain relation although 
orthotropicity is discussed. Only quadratic combinations of displacement 
gradients are retained in the potential energy. A Lagrangian derivation 
is employed with an eigenvalue analysis performed to find the buckling 
loads. A circular or closure element is derived for a static analysis. 
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ELEMENT ID: P“^0 

ELEMENTTYPE: Triangular Plate 



ASSUMED DISPLACEMENT SHAPE: 

Never explicity given. Stated as being: 

Membrane displacement functions are those of constrant 
strain triangle; 

Bending displacement functions are those of the Hsich-Clough- 
Tocher triangle. 


DESCRIPTION: 

This element is an isotropic triangular plate with 3 nodes and 5 DOF 

(u,v,w, B ,Q ) per node. Membrane and bending deformation states 

X y 

are included in the derivation. This is an extension of the HCT triangle 
to postbuckling analyses. Only geometric nonlinearities are considered. 

REFERENCE: [31], [32] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 


R|[PRODUCB^II.n:t u;- 
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ELEMENT ID: P-10 


STRAIN DISPLACEMENT EQUATIONS: 
Strain increment expression 


i uj,: ^ Uk,- ^k,j " j 


u. - displacements in deforrrted equilibrium (initial) configuration 

u. + u. - displacements in deformed equilibrium (incremented) configuration 


DISCUSSION: 

This element employs a linear stress -strain relation with only quadratic 
combinations of displacement gradients in the expression defining the principle virtual dis- 
placements. Convected coordinates are used to form the incremental 
equations of equilibrium and the resulting solutions are related to Lagrangian 
quantities by transformation. The solution procedure is iterative as well as 
incremental . 
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ELEMENT ID: P-11 

ELEMENT TYPE: 

Quadriloterai plate or 
Shallow Shell 



ASSUMED DISPLACEMENT SHAPE: 

Never explicitly given.. Described as 

in -plane: Zienkiewicz-jrons isoparametric quadrilateral [33] 

bending: the Q-19 quadrilateral plate of ref. [34] 


DESCRIPTION: 

This element is an isotropic quadrilateral plate with 4 nodes and 5 DOF 
(u;.v,w, 6^, © ) per node. Membrane and bending deformation states 

are included in the derivation.. This is an extension of the Q-19 plate 
to the geometrically nonlinear regime. 

REFERENCE: [35] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES; 
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ELEMENT ID: P-11 


STRAIN DISPLACEMENT EQUATIONS: (Marguerre Shallow Shell) 

^ 2 \ JiT J " ^ 

<^Vo ^UJ ^uT I f ^tLT\2- J 

^ ^ ^ 2- j T^*- 

yf _ q)' 4' . <j(V 0(UJ (jtJ’ ^ c)^<U/~ 

^‘j 3x 9'A^tj 


DISCUSSION: 

This element employs ci linear stress-strain relation with all combina- 
tions of displacement gradients retained in the strain energy. A 
Lagrangian derivation is employed to obtain the total load-deflection 
equation which is nonsymrnetric and an incremental relation. The non- 
symmetric stiffness matrix does not cause a problem because of the combined 
incremental and iterative solution procedure. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 



ASSUMED DISPLACEMENT SHAPE: 

E THoi W HoJ <q) CJ-;' ^ 

;«/ j-/^- j j j 



(o) _ (o^ 

Hoiu)yoj 


to ti) Ci 3 t /3 (>} 

W,/(x)/Voj Hoi 

vG,^) 



/7^/'}x) = fax^-jixG 

'<^3 = CE"J.<)/j 

DESCRIPTION: 


^ l-lo,lxy--(i-t)/j, y/l 

t^ke^e. OJ^ » -g'y/A , = ^5// J 



This element is an isotropic rectangular plate with 4 nodes and 6 DOF 
(u,v,w, w^, w P®r node. Stretching and bending deformation states 

are accounteJ^ for.^ Transverse deflections are approximated by superposition 
of first order Hermite polynomials and In-plane deflections approximated by 
zeroth order Hermite polynomials. Only geometric non linearities are con- 
REFERENCE: sidered. 

[37] 

VARIATION OF THIS ELEMENT: See also P-2, P-2a, S-5 


ADVANTAGES OR DISADVANTAGES: 
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ELEMENT ID: P-12 

STRAIN DISPLACEMENT EQUATIONS: 



/ 

/^lJ\ 





* 

) 





o>V , 1 

f ' 


/c^^ur) 





/ - Z 




__ c)(X 


c)cj" 




57 ^ 


^ p 

A"-'"" 2? 

<JQ 


DISCUSSION: 

This element employs a linear stress-strain relation with all combinations 
of displacement gradients retained in the strain energy. The third and 
fourth order matrices are evaluated numerically. A Lagrangian derivation 
is employed with a numerical solution with an energy perturbation approach 
used for the post-buckling regime. Initial imperfections are considered. 
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ELEMENT ID: P-13 

ELEMENtTYPE: 

Quadrilateral Plate Composed 
of Four Triangles, i.e., the 
Q-19 Element. 



ASSUMED DISPLACEMENT SHAPE: 

This derivation is carried out in triangular or natural coordinates. The 
displacement function is given explicitly in terms of triangular coordinates 
but the explanation of these coordinates is beyond the scope of this docu- 
ment. The bending displacement function is equivalent to a cubic in the 
rectangular cartesian coordinates. The initial stress stiffness matrix is 
developed for three displacement functions; a compatible cubic, an incom- 
patible cubic, and a linear function. 


DESCRIPTION: 

This element is an anisotropic quadrilateral plate with 4 nodes and 3 DOF 
(w, O ,S ) per node. It is composed of 4 LCCT-11 elements (a version of 

X y 

the HCT element [32]) with the internal DOF reduced out on the elemental 
level. The matrices are derived for a buckling analysis, consequently only 
bending and shearing deformation states are used. 

REFERENCE; [41] 


VARIATION OF THIS ELEMENT; 


ADVANTAGES OR DISADVANTAGES: 

Fully compatible quadrilateral. 
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ELEMENT ID: P-13 

STRAIN DISPLACEMENT EQUATIONS: They are not given explicitly but are equivalent h 



= IT 5 ^ 


DISCUSSION: 

This element employs a linear anisotropic stress-strain relation. Only 
quadratic products of displacement gradients are retained in the strain 
energy in order to develop the initial stress stiffness matrix. A Lagrangian 
derivation is employed with an eigenvalue analysis performed to find the 
buckling loads^ 
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ELEMENT ID: P-14 

ELEMENT TYPE: Triangular Plate 





ASSUMED DISPLACEMENT SHAPE: 

Not given explicitly, but stated as being the BCIZ displacement 
functions of Ref. [43], 


DESCRIPTION: 


This element is an isotropic triangular plate with 3 nodes and 3 DOF (w, <9 , <5^ ) 

X y 

per node. Numerical integration is used in deriving the element, permitting 
thickness variation within an element. The matrices are derived for a buckling 
analysis, thus only bending deformation states are used. 

REFERENCE: [44] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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ELEMENT ID: 

STRAIN DISPLACEMENT EQUATIONS: 
Not explicitly given. 


DISCUSSION; 

This element employs a linear' stress -strain relation in the element derivation. 
Only guadratic products of displacement gradients are retained In the strain 
energy in order to develop the initial stress stiffness matrix. In Ref, 143] 
the stiffness matrix is generated through the use of area coordinates. A 
Lagrangian derivation is employed with an eiganvalue analysis performed 
to find the buckling loads. Ref. [44] presents numerous results using this 
element. 
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ELEMENT ID: P-15 


ELEMENT TYPE: Triangular Plate 


ASSUMED DISPLACEMENT SHAPE: 




DESCRIPTION: 

This element is an isotropic triangular plate with 3 nodes and 6 DOF per 
node. The derivation requires additional DOF of normal slopes at the 
mid-side nodes which are eliminated by imposing a cubic variation in the 
slope. Only bending deformation states are presented. Material non- 
linearities are considered. 


REFERENCE: [48] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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ELEMENT ID: P-15 


STRAIN DISPLACEMENT EQUATIONS: 






Not given explicitly but refers to Ref. [2] for 
the initial stress stiffness contribution 


yj _ y( o ^ 3tcT 

^ ^ x<j i9tj , 5 a 3 ^ 


DISCUSSION: 

This element considers Ufaterjal nonlinearities that obey the Prager-Ziegler 
kinematic hardening theory and therefore accounts for the Bauschinger effect 
for biaxial stress states. The incremental theory of plasticity is used with 
linear elastic and perfectly plastic conditions included as special cases. A 
Lagrangian derivation is employed with only quadratic products of displacement 
gradients retained In the strain energy. Plastic strains are treated in the same 
way as the equivalence between temperature gradients and body forces. This 
results in an ’’Initial strain stiffness matrix" which is not symmetric. A number 
of solution procedures are presented including: 

Predictor Procedure 
Predictor Procedure 
Direct Substitution Procedure 
Direct Substitution Procedure 


Displacement Method 
Strain Method 
Displacement Method 
Strain Methocj 
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Membrane : 

z ) 0-/51 ^ o - v / jh ^ n ./ u , 

3 ) u,/t- fn/i{ v/jL- 0 
H)oA?o 

Vsi‘'$/ex 

O o/_n.j- i/n ^fi = 0 

7) a/j2.-'0 vA " '/k 

8) «A-' - lls/8^ ^.7^.. f s/g;^^ 
DESCRIPTION; 


Bending : 

7) 


s)oj- 


15 

II 

lc>')iLr s -CJ 


'Ooj- 


3>.'^ 


This is an anisotropic parallelogram plate with 4 nodes and 5 DOF 
(UfV,w,(p,6) per node. Stretching and bending deformation states 
are accounted for. Both geometric and material nonlinearities are 
accounted for. The deviation is carried out in natural coordinates. 


REFERENCE: [63] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 


OEBTOA-1^ I’iEi!. - 


nu 
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ELEMENT ID: P-16 

STRAIN DISPLACEMENT EQUATIONS: 

several parts with the membrane and 
form 


A geometrical stiffness matrix is developed for buckling analyses and 
the derivation proceeds from an expression for the generalized load 
matrix. The details of this expression are beyond the scope of this 
handbook. The reader is referred to page 91 of Ref. 63 for the 
derivation of this matrix. 

DISCUSSION; 

This element employs a linear stress-strain relation with only 
quadratic combinations of displacement gradient used to form the 
elemental matrices. A Lagrangian derivation is employed using 
oblique coordinates. An eigenvalue analysis is performed to find 
the buckling loads, and an incremental approach is used in large 
deflection problems. 


The derivation is carried out in 
bending strains having the usual 




n 




Aik 


- ? 


- 2 

inr 
^ . 


- VT* 2 
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ELEMENT ID: P-17 

ELEMENT TYPE: Membrane ^ 

Triangular P3?,te J; 

ASSUMED DISPLACEMENT SHAPE: 



nr Li ^ 5 

DESCRIPTION: 

This element is a plane stress, linear strain geometrically 
ortho tropic triangular plate witlx 6 nodes and 3 DOF (u,v,w) per 
node. Only the stretching deformation state is accounted for. 
Only geometric nonlinearities are considered. 


REFERENCE: 

[65], [64], [63] 

VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGE."!, 
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ELEMENT ID: P-17 

STRAIN DISPLACEMENT EOUATIONS: 



, & i ^ ^ ^ kr 

^ ^5 

DISCUSSION: 

This element has an assumed quadratic displacement state that yields 
a linear strain, linear stress element. A linear stress-strain re- 
lation is assumed and only quadratic terms are retained in the strain 
energy. A Lagrahgian derivation is employed. In Ref. {64] the strain 
displacement relations are modified to include small strain large 
rotation effects 



. / 







^'5 q' ■* 15 ! * 



45 


The derivation includes hybrid 4 and 5 node elements as special cases. 
Centroidal values of the stress resultants are used and the solution 
employs convected coordinates (updated geometry) . The stress-strain 
relation can be orthotropic and nonlinear material behavior can be 
elastic-ideally-plastic , elastic linear hardening or elastic nonlinear 
hardening. In Reference [63] these elements are referred to as the 
TRIM4 through TRIM6 family. 
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ELEMENT ID: p„18 

ELEMENT TYPE: Triangular plate 


ASSUMED DISPLACEMENT SHAPE: 



'Lf = bo bi ?c bi j H 
u_r 'Co Cj X ^ 


4 

- - - 4 C(tj ^ ^ 


— 4 C-a.\ 



DESCRIPTION: 


This is an orthotropic triangular plate with 3 nodes and 12 DOF 

Stretching and 


(u,v.w,e^,8y,©^,€^.€y,e^y,K^,tty.|-<^y) per node. 

bending deformation states are accounted for . 
variation is allowed. 


A linear thickness 


*■ cv,>r/o.,-lor«r 

REFERENCE: [64] 


VARIATION OF THIS ELEMENT: [41] 


ADVANTAGES OR DISADVANTAGES: 

The inplane displacements are complete cubic polynomials and the 
out-of-plane displacement is a complete quintic polynomial. For 
planar structures displacements and rotations will be compatible r 
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ELEMENT ID: P-18 

STRAIN DISPLACEMENT EQUATIONS: 

"‘j ^ 


DISCUSSION: 

This element superimposes the quadratic strain membrane triangle of 
Felippa with quintic displacement bending triangle of Bell. Excess 
degrees of freedom at the centroid - Uj,,V(> are reduced by Gaussian 
elimination and in bending the normal slopes are required to be 
cubic along an edge. The element employs a linear elastic stress- 
strain relation with only quadratic products of displacement gradients 
retained in the strain energy. All stiffness matrix integrations are 
performed using fifth qrder Gauss-Legendre integration formulas. The 
geometric stiffness matrix was developed to be used with an updated 
coordinate system approach. Plasticity is accounted for using Hill’s 
yield criterion and the kinematic hardening proposed by Prager and 
Ziegler. 


= w 1 8 Ux 


A?< 

'3 
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a -J- 
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ELEMENT ID; P-19 

ELEMENT TYPE: dimensional plane 

stress, plane strain, 
or axis3nnnietric ring 
element . 


ASSUMED DISPLACEMENT SHAPE: 

Ui,‘2kk^i j Ua. ^ ^ Ilk^LLa, 

k-i 

where the interpolation functions 
h}^ are given on the next page 



Node r s 

111 
2 -1 1 

3 -1 -1 

4 1-1 

5 0 1 

6-10 

7 0-1 

8 10 


DESCRIPTION: 

This is 3 to 8 node isoparametric element with 2 OOP per node. 

It can be used as an orthotropic plane stress, plane strain or axi- 
symmetric ring element. Only stretching deformation states are in- 
cluded but both geometric and material nonlinearities are accounted 
for. 


REFERENCE; [69] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 




0 





ELEMENT ID 


P-19 


DISPLACEMENT INTERPOLATION FUNCTION 
Define 

R = 1+r 

S = 1+s 

R = 1-r 

s’ = 1“S 

R* = 1-r^ 

S>fc = 1-s^ 


h. = 




-R.S 

iRS* 

i RS* 


Delete if node I is not included 



I - S' 

1*6 

7 

h- = y RS 
1 4 




hj = jls 

-Ci’) H.S- 

■ (■i;)Kc 


h3 = ils 


-(Dk 


^ 



- (i) 


I--S 

■a A, 


-LiA 
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ELEMENT ID: P-19 

STRAIN DISPLACEMENT EQUATIONS: 

Both total Langrangian and updated Lagrangian formulations are 
given on the following page. 


DISCUSSION: 

This is an isoparametric quadrilateral element that can have from 
3 to 8 nodes and thus, model triangular and quadrilateral planforms 
plus refined elements to be used in the transition from coarse to 
fine grids. The element can be used in a linear elastic analysis 
or a nonlinear material analysis. For geometric nonlinearities, 
both total Lagrangian and updated Lagrangian formulations are 
presented . 
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STRAIN DISPLACEMENT EQS. (Continued) 
A. TOTAL LAGRANGIAN FORMULATION 
Incremental Strains 


e 

0 11 

0^1,1 

+ 

0^1,1 

0^1,1 

t 

0^2,1 0^2,1 

+ 

2*- ^0^1,1^ 

+ 

^0^2,1^^^ 

0^22 

0^2,2 

+ 

t 

0’^l,2 

0^1,2 

t 

0'^2,2 0^2,2 

+ 

Ir/ n2 

2*^ ^0^1,2^ 

+ 

^0^2, 2^*^ 

e ■ = 

0 12 

Ir 

2^0“!, 

2 

^ 0“2, 


t 

0^1,1 0^1,2 

+ 

t 

0^2,1 0^2,2 

+ 

t 

0^1,2 O'^l 


t 

0^2,2 

0^2 

.1^ 

Ir 

2 ^0^1,1 

0’^1,2 0^ 

2,1 

0^2,2^ 




0^33 0 ,Q .2 2 \ 0 

( ^1> \ ^l/ 

3u. 3^u. 

it i 

where . = ~Tr~ 5 • = ~7 \ — 

0 1,1 ,.0 0 i,i „0 

8 X. 3 X, 

J J 


(for axis 3 ramietric analysis) 


B. UPDATED LAGRANGIAN FORMULATION 
Incremental Strains 


E- - 
t 11 

t^l,l 

+ 

Ir/ 

+ 

^t^2,l^ ^ 

e 

t 22 

•t^2,2 

+ 

If/ v2 

2^^t^l,2^ 

+ 

^t^2,2^ ^ 

e 

t 12 

Ir 

2^t^l, 

2 

t"^2,l^ + 

1 

2 

i - u^ 

4 1,1 t 1 


H. 


(for axisymmetric analysis) 


where ^u. . = : — 

t i,i „t 

3 X. 

J 
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ELEMENT ID: ^-20 

ELEMENT TYPE: 

Rectangular Plate 



ASSUMED DISPLACEMENT SHAPE: 

The rectangular plate elements are normalized such that x and y range 
from -1 to +1, 

DESCRIPTION: 

This is an isotropic rectangular plate element with 4 nodes and 5 
DOF (u,v,w, w"^',wj) per node. Both geometric and material non- 
linearities are accounted for in the derivation which includes the 
effects of stretching, shearing and bending deformation states. 


REFERENCE: [87], [88] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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ELEMENT ID: P-20 

STRAIN DISPLACEMENT EQUATIONS: 



^ ^13 

DISCUSSION: 

This element was developed for the transient response of a structure 
to blast loads* It uses isotropic materials with the von Mises yield 
condition, the Prandtl-Reuss normal flow rule, and Isotropic work 
hardening. A Gaussian quadrature formula is used for spacial inte- 
gration and control finite differences in the time domain. 



REPRODtJCIBILITr OF THE 
ORIGINAL PAGE IS POOR 


HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P"21 

ELEMENT TYPE: 

"Flat triangular plate-shell element 


ASSUMED DISPLACEMENT SHAPE: 

Transverse displacement field (a) BCIZ, Ref. 143] plate shape functions 
with (b) Razzaque Ref. {91] derivative smoothing [(b) only recently 
added] . 

In plane: linear displacement field; rigid body motions included 

explicity . 


DESCRIPTION: 

This is a triangular plate element with 3 nodes and 6 degrees of free- 
dom/node. Stretching and bending deformation states are included in 
the formulation. 


REFERENCE: [89], [90] 


VARIATION OF THIS ELEM'"i-JT: 


ADVANTAGES OR DISADVANTAGES; 
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ELEMENT ID: p-21 


STRAIN DISPLACEMENT EQUATIONS: 


e 


midplane 

ij 


j_ 

z 


/ <Je-T 

/ 


>C. 


dei~ 





e^. ^ ... strain measured in connected coordinates, a rigid coordinate 

system that rotates with the element. 

rx£ 

u. “ u, where u. and u. are total and rigid body dis- 
i 1 i 1 

placements. 

rotation relative to convected coordinates X, 

Element is applicable to arbitrarily large rotations of plates and 
shells with small strains and moderate variations of rotations within 
an element; a convected coordinate procedure is used. Plate bending 
shape functions are nonconforming and additional incompatibilities are 
introduced by matching of linear in-plane and cubic transverse 
deflections; however, no difficulties have been detected in about 20 
sample problems. 

Second order terms in strain displacement equations are automatically 
omitted whenever variations of rotations within an element is small. 
Element is computationally very efficient and quite accurate because 
convected coordinate formulation implicity includes all rigid body 
modes . 

Element has been used with elastic-plastic (perfectly plastic and linear, 
isotropic work hardening) stress-strain law in a transient analysis 
program with explicity temporal integration and (unreported) a static 
program. 


u . 

1 


def 


ki ' ' ’ 

DISCUSSION: 
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ELEMENT ID: P"22 


ELEMENTTYPE: Triangular membrane 


ASSUMED DISPLACEMENT SHAPE: 
Linear 



DESCRIPTION: 

No. of nodes : 3 

Degree of freedom: 9 

Coordinate system: Cartesian 

Isotropic materials 


REFERENCE: [114], [113] , [107] , [115], [117], [118] 


VARIATION OF THIS ELEMENT: Jl] , [2] , J21] , [60] , [94] 

ADVANTAGES OR DISADVANTAGES: 
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Advantages or disadvantages : 
Advantages: Computational simplicity 


Strain Displacement Equations ; 




+ u . . 
D /I 


+ u .u . ) 


Discussion : 

a. Stress-Strain equations: ■ + y— ) + pG^^ 

b. Terms retained in the strain energy: All 

c. Derivations from other than strain energy considerations: 

Galerkin Methods, Conservation of Energy * 

d. Unique derivation techniques: See (c) above. 

e. Lagrangian, Eulerian, convected: Lagrangian (convected on 

material) 

f. Numerical integration schemes: 

For static problems: Incremental Loading (Euler's Method) 

plus Newton-Raphson corrections. 

g. Required or preferred solution procedure: 

For dynamic problems: Incremental Loading. Divided central 

difference methods plus a Lax-Wendroff method. 
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ELEMENT ID: P-23 



DESCRIPTION: 


No. of nodes : 
Degrees of freedom: 
Coordinate system : 
Isotropic materials 


6 

18 

Cartesian 


REFERENCE: » IlUX, IIO 7 J, Jl 1 5] , IH7J , IUBJ. 


VARIATION OF THIS ELEMENT: I63J , I64J , IBSJ . 


ADVANTAGES OR DISADVANTAGES: Advantages: Computational slmpltcrty 
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strain Displacement Equations ; 


Y.. = 7t(u. . +u. . + u .u .) 

13 2 j ,1 m,i m,3 


Discussion ; 

a. Stress-Strain equations: ^ ^ 

^ij ®^ji 

b. Terms retained in the strain energy; All 

c. Derivations from other than strain energy considerations; 

Galerkin Methods , Conservation of Energy ' 

d. Unique derivation techniques; See (c) above. 

e. Lagrangian, Euler ian, convected: Lagrangian (convected on 

material) 

f. Numerical integration schemes: 

For static problems: Incremental loading (Euler's Method) 

plus Newton-Raphson corrections, 

g. Required or preferred solution procedure; 

For dynamic problems: Incremental Loading. Divided central 

difference methods plus a Lax-Wendroff method. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 
ELEMENT ID: P“24 


ELEMENT TYPE: TRIM3 

(triangular membrane element in 3 space) 



Pi. 


DESCRIPTION: Number of nodes : 3 

Degrees of freedom: u, v, w at each node, all together 9 

Deformation type : stretching, shearing 

Isotropic or anisotropic 


REFERENCE: 


♦ 

. VARIATION OF THiS ELEMENT: 111, l2j, l2lj, I60Jv Il07J, IH3J, 111^1,1113,11171,1)18]. 

For application in large displacement dnd large strain the special'element TIIIMP3 
was derived from TRIM3- This applies also to higher order elements of the TRIM 
family, which are not listed here (as TRIM^ and TRIMP6, TRiMlO and TRIMPIO). 
ADVANTAGES OR DISADVANTAGES: 

D I SCUSS ION: Based ‘on the concept of natural stress- invariants for large rotations. 

Normally the application rate is low. But considering problems of large straijis 
and large displacements this element is superior to even higher order element's. 

Ref. I95].' 


108 


REPRODUCIBILITY OF Til 
ORIGINAL PAGE IS POOR 


HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: P"25 



Flgure i 






IH =0 
\ 


t/2 


t 








' ' 



DESCRIPTION: Number of nodes : 3 

Degrees of freedom: u, v, w,' ^ ^ , <j> , at each node 

X y ^ 

(<j). Angle of rotation), all together 18 
Deformation type : stretching, shearing, bending without shear 

Isotropic or Anisotropic 


REFERENCE: , {38J 


VARIATION OF THIS ELEMENT: An additional centrifugal stiffness part Was 

developed for the application to rotating parts. 

ADVANTAGES OR DISADVANTAGES: Advantages ; Simple in application and use. No 
numerical integration, very fast element processing. Good experience in buckling 
and large displacements. 

Disadvantages ; Not fully compatible for non-plane 
structures. Curved structures have to be- approximated by faceting thus involv- 
ing a relatively large number of elements. 
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ELEMENT ID: P-26 


4 /v/ 






/ 


‘B,ux 


DESCRIPTION: Number of nodes : 

Degrees of freedom : 


REFERENCE: [99]. 


Deformation type: 


Vertices 

w w, X w, y w, XX w, xy w, yy 
3 4 5 6 7 8 

Midside nodes: 
w, n • ' 

9 

t 

n stands for a local direction normal to the 
edge;, number beneath the symbols indicate 
the internal number of the corresponding degree 
of freedom; all together 21 
bending 


Isotropic or Anisotropic 

VARIATION OF THIS ELEMENT: Elimination of midside nodes leads to TUBA3. 


ADVANTAGES OR DISADVANTAGES; 

Advantages : More accurate then TRIB3 (based on the same effort) , fully 

compatible. 

Disadvantages: Jumps in thickness provide difficulties. _ 

DISCUSSION; Based on the concept of natural stress invariants for large rotations. 
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ELEMENT ID: 

ELEMENT TYPE: 

(See next page also) 

ASSUMED DISPLACEMENT SHAPE: 

LLf - CL, + CXz f 

^ 2 . " ^3 f ^ f ^ y 

Geometry: 

?" ^ O' f }, o., = f 

DESCRIPTION; 

This element is a doubly curved, isotropic shell of revolution restricted to 
axisymmetric loads. It has 2 ring nodes and 3 DOF (u^,U 2 ,^) per node. 

Membrane and bending deformation states are included in the derivation. 
Geometric and material nonlinearrties are considered. 

REFERENCE: [42] 

VARIATION OF THIS ELElVjENT; [46], [47] 

ADVANTAGES OR DISADVANTAGES; 
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BEPRODUCIBILITY 
ORIGINAL PAGE I 


ELEMENT ID: S-1 

STRAIN DISPLACEMENT EQUATIONS: 




6, - e. + 


£ . 

2. ) 




•5 




eR< 


u .. V f ^ 

^ o ■A. 



e< ~u 


tj 


JS 


X. 


r '5 


es> ^ 


(uC-OS ^ <j> ) 






DISCUSSION: 

This element considers material nonlinearities which obey the won Mises 
yield criterion and isotropic hardening law. The incremental theory of 
plasticity is used. The linearly elastic and perfectly plastic conditions are 
included as special cases. The mechanical sublayer model is used for the 
elastic -plastic analyses. A Lagrangian derivation is employed and only 
quadratic products of displacements and/or displacement gradients are retained 
in the virtual work expression. The solution is obtained by an incremental 
stepwise procedure. 
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ELEMENT ID: 


S-la 


ELEMENT TYPE: Doubly Curved 

Sbel! of Revolution 






/.a 



REFERENCE: [46], [47] 

VARIATIONS FROM BASIC ELEMENT: S -1 

This element uses the same displacement function as element S~l, but 
the following form of the strain -displacement equations is employed 


cJc, 

'"a 


-h 




iU _ ^ _ S 


Uc^C-k ) , ^,7 / ' 5 ^^^ U.C. \ 

C? S ^ ‘Z )T ^ f 


u 


Ue. 

A 


Material noniinearities are accounted for by the incremental theory of 
plasticity and the von Mises yield criterion. It is assumed to work harden 
according to an isotropic strain hardening criterion. The mechanical sublayer 
model is employed. Numerical integration , is used to form the stiffness matrices 
A direct solution using a sub-matrix form of the Crout-Banachiewicz reduction 
procedure was developed. A Lagrangian derivation is employed with all 
products of displacements and/or displacement gradients retained in the virtual 
work. 
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ELEMENT ID: S-2 

ELEMENT TYPE: Doubly curved 

Shell of Revolution 


ASSUMED DISPLACEMENT SHAPE: 


a," o(, f o(j 

Uz= 

<5l = f'’ 

^ = °</ 2 . ^ 

DESCRIPTION: 


7 = 



Bottom Top 



This element is a doubly curved, isotropic, three layer sandwich shell of 
revolution. It has two external ring nodes with 5 DOF (u^,U 2 , 3^) 

per node and three internal ring nodes with DOF (u^, ( and (u^ 

respectively. Membrane, bending, and shear stiffnesses are included in the 
derivation. Geometric and material non linearities are considered. 


REFERENCE: [49] 


VARIATION OF THIS ELEMENT: S-2a 


ADVANTAGES OR DISADVANTAGES: 

It may be difficult to incorporate this element in a general purpose program 
because of the non-typical nodal degrees of fnsedom. This element satisfies 
completeness and compatibility requirements. 
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REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


ELEMENT ID; 


S-2 


STRAIN DISPLACEMENT EQUATIONS; 


f" ^ -- fe" ? ^ ^ .1 

Co ^ Jk 

•‘•rc v'f, Ki^ 3 

rf^Kzc'frf. /Li/J 



H_= / 5 ;j 

i ~- 



IK /Lie, 

K-la .Ktc. 





-+A' 




r (<^/c " 2:(^)x ^ 6c - ) 

g* ei - 2 ( Tt ^ Kif^^zc - ^ [^c. 6c^ 0: ^ ^ i 9 

U,^^Kje,c ^ 

C ^ZH, ^ Kz.<i) ^Z<L *" Zr S'lVl^^ 


c *Ku/| ^ “ 2t^ (X” ^ 


6>J loaves. » ^ 


VC 


'I C«J 


o^U, 

<?5 ’ R, J 

()cJ 




6. 


(j> ^arsii'] !),k 


^ > ^.a "--^ 

H 


/< 


0^ 


2C 


dC’5 


K 

N ) ' '•W 

A =-^ 

I >^\4 J5 


2^. . 

r (Lc?5 








X= ^ ^ > K 

'Ol- ■^ ) ' ^ ) r - Rxd^ 

This element considers material non linearities which obey the von Mises yield 
criterion and isotropic hardening lav/. The incremental theory of plasticity is used. 
The core material is elastic whereas the facings may be inelastic. An incremental 
variational form based on a moving reference configuration is used to derive element 
equilibrium equations. The incremental Lagrangian forms of strain displacement 
equations include terms corresponding to large rotations. Only quadratic products 
of displacements and/or displacement gradients are retained in the virtual work. 

The solution is obtained by an incremental stepwise procedure with post-buckling 
solutions obtained by augmenting the structure to obtain a positive definite stiffness 
matrix* 
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ELEMENT ID: S-2a 




ELEMENT TYPE; Doubly curved 
shell of revolution 
cop element 


REFERENCE: {49] 



VARIATIONS FROM BASIC ELEMENT: S-2 


This is the frustum element of S-2 specialized to a cap element. The 
assumed displacement shapes have the form 


jj-h 
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ELEMENT ID: 


ELEMENT TYPE 


ASSUMED DISPLACEMENT SHAPE: 

Iff ; • ^ < 

^ ^ <loS I S ^ ( '^5- ^ ^ Si I © 

;--o 

^ ^ 0 (q s)si^i e C^y s)(i.05t© 

; = o * 

1)3 _ 

S^)(Lasi © ^ ^ ('C<j 

I TO i‘0 

DESCRIPTION: 

This element is a doubly curved orthotropic shell of revolution developed for 
arbitrary nonaxisymmetic loads, it has 2 ring nodes and 4 DOF (u,v,v/, rotation) 
at each node. Membrane and bending deformation states are included in the 
derivation. Only geometric nonlinearities are considered. 


REFERENCE: [50] 


'^e<^5 t sir 


( • * I 

- ^ C f 0(-^S + 0(J -I- o(^ 



VARIATION OF THIS ELEMENT: S-3a, S-3b, S-3c 


ADVANTAGES OR DISADVANTAGES: 

In a later paper [53] .the authors show that fourth order products of displace- 
ment and/or displacement gradients should be retained in the strain energy to 
improve solutions. 
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ELEMENT ID: 


S-3 


str'ain Displacement equations: 

I A Z. A 


/I 2. 


^ ^ , I AZ. A I 'I 


V 

ds J 

c)u 


i_ <^e-z 3 

~r 


& 

/ A 


^ 5in ^ i 


/I A 

>3 ^2i 

\ <^e 
^ 9& 


+-^ 


r 5/ 




A 

■~gj 


I f j_ r , / / “1 

^5 ” ^ “ r- ^ “ r L ^ USmf i-CJ(lo^<p ] 

4 I c)iK V . / c?l/ A i5cj“ 


^ j , A t_ cP«r 

As ) ^ 13 ' j ^2s'‘ y~ 


ii 




(Lo^ 


/ 


DISCUSSION; 


This element employs a linear stress-strain relation in the element derivation. 
Cubic products of displacements and/or displacement gradients are retained 
In the strain energy, but are treated as an additional force column. The 
element is developed through the use of the non-shallow shell theory of 
Novozilov. The displacement function is that of [51] as extended in [52]. 
The stiffness matrix is evaluated through numerical integration over the 
length and the arbitrary loading represented by a Fourier series in the 
circumferential direction. A Lagrangian derivation is employed with the re- 
sulting equilibrium equations solved by one of three methods: load increment, 
iteration, or a combination of the two. 
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ELEMENT ID: S-3a 

ELEMENT TYPE: Doubly curved 

Shell of Revolution 



5-3) 


REFERENCE: [53] 

VARIATIONS FROM BASIC ELEMENT: S-3 

This element is identical to that of S-3, except that the fourth order 
products of displacements and/or displacement gradients are retained in 
the strain energy. Several examples are included that illustrate the 
improved solutions. 
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ELEMENT ID: 


S 


ELEMENT TYPE: Doubly Curved 

Shell of Revolution 



REFERENCE: [61] 

VARIATIONS FROM BASIC ELEMENT: S-3 

This element apparently uses the same cubic displacement function as S-3 
but the generalized displacements are given explicitly as 

U- ^ t) I Z~\ 

"V ~ X d ^<55/(9 

/ = <? ta f 

The w displacement is assumed to be 

aj“|I T J 

in the strain energy contributions introduced by nonlinear theory, and 

The element Is a doubly curved orthotropic shell of revolution developed for 
arbitrary nonaxisymmetric dynamic loads. In other respects, it is similar to 
element S-3 except that fourth order products of displacements and/or displacement 
gradients are retained in the strain energy as recommended in element S-3a. 
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ELEMENT ID: S~3c 

ELEMENT TYPE: Doubly Curved 

Shell of Revolution 


j 1 

1 
1 

REFERENCE; [60] 

VARIATIONS FROM BASIC ELEMENT: S-3 

This element uses the cosine harmonic displacement function as shown in S-3 

If? 

(i.e., H terms), but obtains an eiganvalue solution for the buckling load. 

It is a doubly curved isotropic shell of revolution developed for nonaxisymmetric 
eiganvectors. It has 2 ring nodes and 4 DOF (u,v,w, rotation) at each node. 

Membrane and bending deformation states are included in the derivation. Only 
geometric nonlinearities are considered. The Sanders' strain displacement equa- 
tions for; moderate bending are given as 

6 ^^ - (^17) (t'Vs f E % 

where % - - ( ^ ^ '‘(r) (j>) ^ ^ ) 

By dropping the terms with asterisks, the Donnel I -Mushtari -Vlasov strain- 
displacement relations are obtained. A linear stress-strain relation is used. 

A Lagrangian derivation is employed with quadratic products of displacements 
and/or displacement grqdiepts retained in the strain energy. 
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ELEMENT ID: 


S-4 


ELEMENT TYPE: Conical Shell 

of Revolution 



ASSUMED DISPLACEMENT SHAPE; 

For convected coordinates, the deforrrration displacements are the 
nodal rotations relative to the x axis 






^z. 

A 

+ a, X 


DESCRIPTION; 

This element is an isotropic conical shell element developed for a large 

displacement small strain, elastic-plastic dynamic problem. It has 2 ring 

nodes and 3 DOF (u ,u ,(p) at each node. Convected coordinates are 
X y ' 

used to derive the elemental matrices. Membrane and bending deformation 
states are included in the derivation. Both geometric and material non- 
linearities are included, 

REFERENCE: [54] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 

The method is stated as being computationally efficient. 
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element ID: S-4 

STRAIN DISPLACEMENT EQUATIONS: 

A «/«/' , ^ /A\ 

r ^ f ^ >o \ 


DISCUSSION: 

This somewhat brief paper indicates the element can account for material 
non linearities by using numerical quadrature through the thickness. Be- 
cause convected coordinates are used, the geometric non linearities that 
arise are accounted for by transformation between the convected and global 
systems. The strains are linearly related to deformation displacements 
relative to the convected coordinate system. The dynamic equations of 
equilibrium are solved by a central difference explicit integration procedure 
which uses total nodal forces. 
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ELEMENT ID: 


ELEMENT TYPE: 


S-5 

Cylmdrical 
Shell Elemenf of 
Rectangular Plan- 
form 





ASSUMED DISPLACEMENT SHAPE: 

2 . a 


cj(^ re)-T ^ \ ^ -+ 

/r/ )';/ •- */ -i ' V ^ 

dl) V s “1 

CU-. //"( 5 )= 

Similarly for the O direction, replace s by r <9 , and I by b. Similar 
expressions for u and v. 

DESCRIPTION: 

This element is an isotropic cylindrical shell element of rectangular planform 

with 4 nodes and 12 DOF (u,v,w,u ,v ,w ,u ,v ,w ,u ,v ,w ) per node. 

XXX y y y xy xy xy 

Membrane and bending deformation states are accounted for. Displacements 
are approximated by products of one-dimensional cubic Hermite interpolation 
polynomials. Only geometric nonlinearities are considered. 


REFERENCE: 


16] 


VARIATION OF THIS ELEMENT: See also P-2, P-2a, P-12, S-5a 


ADVANTAGES OR DISADVANTAGES: 

The bicubic interpolation functions for u,v, and w admit to six linearly 
independent displacement states of very little strain energy. Consequently, 
this element resolves the rigid body mode problem. The selected nodal DOF 
may not be compatible with other elements in a general purpose structural 
analysis computer program. 
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ELEMENT ID: S-5 

STRAIN DISPLACEMENT EQUATIONS: 


2 ( 7TJ z ( -jpr/ 

^ - _L_ tL 4- _L_ / ^ ^ 

® r* <90 r* 2.r^\^&/ r^. ( 

55 r 5<9 r 5 j jr ( ^5 <50 <-75 / 


DISCUSSION: 

This element employs a linear stress -strain relation with all combinations 
of displacements and/or displacement gradients retained in the strain 
energy. A Lagrangian derivation is employed with a numerical solution 
obtained by direct minimisation of the total potential energy. 
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ELEMENT ID: S-5a 

ELEMENT TYPE: Doubly curved 

Rectangular Shell 

REFERENCE: [55] 

VARIATIONS FROM BASIC ELEMENT: S -5 

This element employs the displacement function of Ref. [16] for a 
doubly curved shell element. Only geometric nonlinearities are considered. 

The shell strain displacement equations are 

Y'.uy ^ ^ (r,U^ ) 

+ a' r^c^)(cj^) - Z (rg 

where r^ == 1/R^ and r^ = constant. The curvilinear coordinates 

are ^ with subscripts denoting derivatives. A systematic integration 

method is developed that simplifies the element formulation. It can account 
for a non-uniform distribution of initial stress. Apparently cubic products of 
displacements and/or displacement gradients are retained in the strain energy, 
but since an eiganvalue solution is performed to determine dynamic modes and 
frequencies a stress state is assumed. As used in the context of this document, 
this is equivalent to an initial stress state with quadratic product of displace- 
ments and/or displacement gradients retained in the strain energy. A Lagrangian 
derivation is employed with an eigansolution performed to find dynamic modes 
and frequencies. The linear stiffness matrix is derived in Ref. [56] based on 
the strain displacement equations of [57]. 



126 


HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: S-6 

ELEMENT TYPE: Doubly Curved 

Rectangular 
Shell 


ASSUMED DISPLACEMENT SHAPE: 

U-- (^0^, (0(2,“^ } ^/ ■" ^ 6(2 ^2. Uj ■■ “Cij ( 5 < 2. U«yJ 

V- (o^, \/^ - ~ (oi.^-C^) ] /o.;4 

oS ~ the product of one dimensional cubic Hermite interpolation polynomials 
similar to that of 5-5. It is given explicitly on the following pages, 

DESCRIPTION; 

This element is an isotropic doubly curved shell element of rectangular 
shape with 4 nodes and 6 DOF (u,v, w, ^o<^, (^ 2 ^ 0 C 2 ) node. 

Membrane and bending deformation states are accounted for. Displacements 
are approximated by Lagrangian and Hermitian interpolation polynomials. 
Only geometric nonlinearities are considered. 

REFERENCE: [58] 


VARIATION OF THIS ELEMENT: PhS 


ADVANTAGES OR DISADVANTAGES: 

This element violates both the postulated requirements that rigid body motion 
be included in the assumed displacement shape and that displacements be con- 
tinuous across interelement boundaries. The authors found these factors had 
little influence on the solutions obtained. 
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ELEMENT ID: S-6 

STRAIN DISPLACEMENT EQUATIONS: Extension of Novozhilov equations [54] 

Deep shell eqs. given on following pages. 


DISCUSSION: 

This element employs a linear stress -strain relation with only quadratic 
products of displacements and/or displacement gradients retained in the 
strain energy. The element is developed through the use of an extension 
of the shell equations of Novoshilov [59], Explicit forms of the stiffness 
equations are obtained in related references. A Lagrangian derivation is 
employed with an eigansolution of the resulting equations. A Guyan type 
reduction is applied to the initial stress stiffness matrix on the structure 
level and the resuits compared to unreduced solutions. 
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ELEMENT ID: S-7 

STRAIN DISPLACEMENT EQUATIONS: 

'(Musbtar i ' s Eqs . ) 

DISCUSSION; 

This element employs a linear stress-strain relation with all combina- 
tions of displacement and displacement gradients retained in the strain 
energy. A Lagrangian derivation is employed with a numerical solution 
obtained by a perturbation technique. The neutral equilibrium state i.s 
considered as is the solution through a critical point. 
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ELEMENT ID: S-7 

ELEMENT TYPE: 

Doubly curved .shallow 
triangular shell 



ASSUMED DISPLACEMENT SHAPE: 

w = quintic polynomial in dimensionless triangular coordinates 
u = cubic polynomial 
V = cubic polynomial 


DESCRIPTION: 

This element is an orthotropic doubly cuTrved shell element of triangular 
shape with 3 corner nodes and 12 DOF (u,u,x,u>y,v,v,x,v,y,w, w,xjW,y 
w,yy) . A centroidal node with 2 DOF (u,v) is eliminated by static 
condensation. The normal slope is constrained to a cubic variation 
along the edge. Membrane and bending deformation states are accounted 
for. Only geometric nonlinearities are considered. 

REFERENCE: 162] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 
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ELEMENT ID: S-8 


ELEMENT TYPE: Doubly curved 

Shell of Revolution 


ASSUMED DISPLACEMENT SHAPE: 

, , \ . I II _ I 


CO, 



'■IT ■ ■(- f^oi < 

where Hot^Cf), HoiCO, are cubic Hermitian interpolation 

polynomials 

• Geometry ; 1 |_t ( 0. f 3 f % 0 ^ C3 f 2 A - 2 ^ V 4 C f f V 

DESCRIPTION: 


This element is a doubly curved, orthotropic shell of revolution 
restricted to axisymmetric loads. It has too ring nodes and 6 DOF 
(u,w,'X,e°, 

V per node. Membrane and bending deformation states are included 

in the isoparametric derivation. Geometric and material nonlinearities 
are considered. 

REFEREfVCE: [64] 


VARIATION OF THIS ELEMENT: [67] 


advantages or DISADVANTAGES: 


f: 
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ELEMENT ID: 3-8 


STRAIN DISPLACEMENT EQUATIONS: Sanders theory for small strains and 

moderate relations. 

€-e= ? ka 

^ ? kse 

3i ^ f, ' < ^ -r ^ - If , k, = -^ y 




- 1C 

Rt, ^ 


DISCUSSION: 


This element employs a linear elastic stress-strain relation with only 
qudratic products of displacement gradients retained in the strain 
energy. The geometric stiffness matrix was developed to be used with 
an updated coordinate system approach. Integrations are performed 
numerically with up to eight order integration formulas available. 
Plastic material properties obey Hill’s yield criterion. Up to 20 
layers through the thickness may be used to calculate effective plastic 
forces and moments. 
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ELEMENT ID: S-9 

ELEMENT TYPE: Axisymmetric Shell 


ASSUMED DISPLACEMENT SHAPE: 



A displacement shape is not assumed. The exact stiffness matrix, 
including nonlinear effects, is computed by using Fourier series 
expansions in the circumferential direction and a numerical inte- 
gration procedure to solve the resulting differential equations. 


DESCRIPTION: 

Each element has two circumferential line nodes with four degrees 

of freedom per harmonic at each node, u , v , w and : 

^ n n n 0n 

(Continued next pg.) 


REFERENCE: [73], [74], [75], [76], [77], [78], [79]. 


VARIATION OF THIS ELEMENT: 

An element in the form of a circumferential segment of a torus with 
meridional lines nodes is used in Ref, [80]. 

ADVANTAGES OR DISADVANTAGES: 

The advantages of this formulation are many. First, the use of the 
numerical integration procedure, instead of an assumed displacement shape, 
yields "exact" stiffness matrices. Then, the transformation of the influence 
coefficient matrices into stiffness matrices, which is unique (other formula- 
tions using numerical integration perform all computations with the influence 
coefficients), allows the use of methods developed for finite element technology. 
These are; arbitrary branching, use of equilibrium corrections, use of Guyan 
reduction scheme and the method for solving eigenvalue problems. 
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ELEMENT ID : S-9 


DESCRIPTION : 


This program can analyze orthotropic thin shells of revolution > 
subjected to unsymmetric distributed loading or concentrated line loads, 
as well as thermal strains. Furthermore, a shell with arbitrary boundary 
conditions, under loads which vary arbitrarily with position and under a 
temperature variation through the thickness, is tractable with this program. 
The shell can consist of any combination of the following geometric shapes; 


1 . 


2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 


Ellipsoidal - spherical (offset from the axis of revolution 
allowed) 

Ogival - toroidal 
Modified ellipse shape 
Conical - circular plate 
Cylindrical 

General point input geometry 

Dummy geontetry slot to be filled by the user 

Discrete pin^ 

Elastic support . 


The shell wall cross-section can be a sheet, sandwich, or reinforced sheet or 
sandwich. The reinforcement can cons. J of rings and/or stringers, a waffle 
construction rotated at any angle to the principal coordinates, or an isogrid 
construction. General stiffness input options are also available. The rein- 
forcement material properties can differ from those of the main shell, and a 
temperature variation can cause different properties in the two face sheets 
of a sandwich shell . 
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ELEMENT ID: S-9 

DESCRIPTION: 

The basic approach to the problem is to cut the structure into 
several shell regions. These regions need to be singly-connected shells, 
and can only have line loads applied at their end points. There are no 
restrictions on geometry, or uniform or thermal loads. The regions are 
further subdivided into several shell segments, each being free to have its 
own geometric shape, provided that the shape falls into one of the categories 
mentioned above. 

There is a considerable latitude in what can be done within each shell 
segment. The thickness of any segment can be s 3 nnmetrically tapered and it 
can contain up to 14 points of discontinuity, provided that the segment center- 
line remains continuous and describable by a single shell geometry. A tem- 
perature distribution through the thickness can be specified at three points 
in a homogeneous shell, and 4 points in a shell of rigid core sandwich con- 
struction. The distribution is considered to be linear between these points. 
Thus, it is possible to approximate temperature distributions other than 
linear distributions. In the event of physically discontinuous shell center- 
lines, a kinematic link is available for use in the analysis. The link 
relates displacements across the discontinuity. This link may be used between 
regions, and between segments within a region. Discrete offset rings are 
also available for use within or between regions. 
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HANDBOOK OF NONLINEAR FINITE ELEMENTS 
(Continued) 


ELEMENT ID : S-9 


DISCUSSION: 

This element is available in the STARS (Shell Theory Automated 
for Rotational Structures) computer programs developed at Grumman Aero- 
space Corporation and partially supported by NASA/MSFC. The following 
STARS programs are available: 


1. 

STARS 2S 

(Static Analysis) 

2. 

STARS 2V 

(Vibrations) 

3. 

STARS 2B 

(Buckling) 

4. 

STARS 2P 

(Geometric and Material Nonlinearities) 


The use of an accurate shell theory to analyze structural shell 
problems usually involves complex mathematics and numerical techniques, 
which are nearly impossible to treat without the aid of automated procedures 
On this basis, a digital computer program (STARS) based upon the Love- 
Reissner first order shell theory has been developed. 

The shell equations are first cast (after Fourier expansion in the 
0 coordinate) in the form of eight coupled first-order equations of the form 

[^•■“^7)], y = (1) 

where (/?, 6 are the meridional and hoop coordinates, respectively, and n 
signifies the Fourier harmonic. These equations are then integrated by 
using some standard procedure such as the Runge-Kutta method. Eight in- 
fluence coefficient solutions are thus obtained corresponding to eight 
independent initial condition vectors, and a particular solution for zero 
initial variables may be obtained for each independent loading pattern to 
be investigated. Due to the numerical difficulties caused by exponential 
growth or decay of the influence coefficient solutions, the shell must be 
segmented into pieces of limited size to obtain a satisfactory solution. 

Thus for each shell segment a set of (8+P) solutions is obtained, where F 
is the number of individual load patterns being considered. The total sets 
are then combined to satisfy the segment continuity conditions, and the 
overall boundary conditions for the shell, thus yielding P solutions for 
stresses and displacements in the case of a static analysis. 

These solutions represent influence coefficients for the shell seg- 
ment which are then transformed into stiffness and load matrices. Finite 
element procedures used in the direct stiffness method are then applied to 
obtain a solution. 
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(Continued) 


ELEMENT ID : S-9 


DISCUSSION: 

For the solution of buckling problems a new approach is taken. It 
is first recognized that the stability or vibrations problem for the shell 
is actually transcendental in the eigenvalue. Then, in successive passes 
through the program, the elastic stiffness matrix of the shell (K^) , and the 

prestressed shell stiffness matrix [K^CX)], (which includes the nonlinear 

terms) for an assumed prestress value, are calculated from the differential 
equations. The two matrices are subtracted to isolate a matrix containing 
the effect of the prestress: 

[Kj(A)j = [KpU)] - [KqJ (2) 

and a linear eigenvalue approximation is formulated: 

^i 

([Kq] [K^(U-1)J)[A] = [0] (3) 

i-1 

where [A] is the deflection vector of the segment ends. As evident from the 
notation, Eq . (3) is solved by a series of iterations, with convergence 
being signaled by this fashion the nonlinearity of the eigen- 

value problem is retained. The iteration procedure has been found to exhibit 
quick convergence. A single iteration is equivalent to the process utilized 
in most finite element programs, whereas further passes include the refine- 
ments of accounting for the nonlinear dependence of the incremental stiffness 
matrix upon the eigenvalue. The formulation of Eq. (3) includes the effects 
of all 'consistent' nonlinear terms including predeformation. 

An incremental Lagrangian approach is used for the geometric nonlinear 
analysis. The nonlinear shell formulation is based on a Love-Reissner-r ' 
Kempner shell theory and is valid for small strains and moderate rotations. 
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ELEMENT ID: S-10 

element TYPE: 

Doubly curved 
shell of revolution 



o-*i node CQ 


ASSUMED DISPLACEMENT SHAPE: 



* 1 i H; 


a I 


-Slv^0; 
, >^3 <=)t 


oCu 


£M\d 

r 
. ? 





=(>=,3CfJ = ^0+3f)0-f’) 

DESCRIPTION: 

This is a doubly curved, isotropic, shell of revolution. It has four 
nodes, two internal, with 3 DOF (u,w,a) per node. This is the 
degenerate isoparametric shell element of Ahmed, Ref. [83]. Membrane 
and bending stiffnesses are included in the derivation. The partic- 
ular application includes geometric nonlinearities in a viscoelastic 
shell . 


REFERENCE: [82] 


VARIATION OF THIS ELEMENT: In Ref, [84] this element is used in a nonlinear 

elastic-viscoplastic analysis, 

ADVANTAGES OR DISADVANTAGES: 


RliJPRODUCIBiLlTY OF OM 
ORIGINAL PAG® P 
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ELEMENT ID: s-10 


STRAIN DISPLACEMENT EQUATIONS: 




where the superscript 1 refers to 
current configuration 


DISCUSSION: 


the 


This element was developed by generalizing the infinitesimal theory 
of linear viscoelasticity to small strain large relation applications . 
The form given in this application is in terms of the Lagrangian strain 
rate and the second Piola-Kirchhoff stress tension. The viscoelastic 
material is characterized using Preny series expansion for the rela- 
tion modulus. A Lagrangian strain description is used with both con- 
servative and nonconservative loads derived using the virtual work. 

The resulting equilibrium equations consist of the incremental and 
geometric stiffnesses. 
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X 3 

a- CL, 4 4 cLyy)T^4o.< -.arjf 4 0-g'^^rj^ +o^'^r^^4o.,Qix^ 

'XT- <^1> 4 <3.,7.t -^O-K^ ■''®i 5 f 40-1, t. 4 0 -vT) f ^ < ^'8 f 4 0.>c\^' ^ 

C 0-j^^ 4.CK^^.■^ -VO-x^r^^Ci-ai^^ HG-aT)f 4<V-2.Bf *7. 4 aa.^^ 0-jo 

DESCRIPTION: 


This element is an isotropic doubly curved, shallow, triangular shell 

element. It has three vertex nodes with 12 DO? (u,v,w,u ,v ,w ,u v^,w^, 

w ,w ,w )per node and a controidal node with two DOF that 

^x]’ 

are removed by static condensation. It is developed in a cartesian 
coordinate system and includes both bending and stretching deformation 
states. 


REFERENCE: |86.1 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES; 

The element uses high precision displacement functions resulting in 
better accuracy with fewer elements than for elements with lower order 
displacement functions. 
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ELEMENT ID: S-11 

STRAIN DISPLACEMENT EQUATIONS: 
middle surface strains 

~ ^ 2 . 

bending strains 
DISCUSSION: 

This is the Cowper element of Ref. [87], extended to the large 
deflection regime. Linear elastic stress strain equations are used 
and all terms up to fourth order are retained in the strain energy. 

The formulation being based on the above strain-displacement equations 
is a Lagrangian formulation. The integration is performed exactly. 

The Newton-Rophson procedure is used to solve the equations of 
equilibrium. 
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ASSUMED DISPLACEMENT SHAPE; Cubic in area coordinates for 
membrane and bending displacements. 


DESCRIPTION: Uses orthogonal curvilinear coordinate system to describe 

doubly curved geometry of element, with user-input coordinates and curvatures 
at the three corner nodes . Degrees of freedom are three displacements and their 
derivatives with respect to the curvilinear coordinates at the three corner nodes 
and the displacements at a central node, a total of 30. Constraints are imposed 
to obtain approximate slope continuity, ‘using a Lagrange Multiplier Technique. 
Retains the Lagrange Multipliers explicitly, in the set of unknowns to be solved. 

REFERENCE; IH9] . 


VARIATION OF THIS ELEMENT; 


ADVANTAGES OR DISADVANTAGES: Advantages include simplicity due to cubic dis- 

placement forms, with competence due to using cubics for the membrane displace- 
ments . ' 

Disadvantages include the need to handle Lagrange Multiplier and the fact that# 
rigid body motion causes small element straining. The element shows accuracy' 
in applications. 
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ELEMENT ID: S-12 


STRAIN displacement EQUATIONS: Green’s strain in a La^rangian formulation 
is used. Effectively, the nonlinear terms correspond to the flat membrart: 
equation 


8u . 1 
■" 3x 2 




etc, , 


the membrane nonlinear term being dropped. Actually, a deep shell strain- 
displacement formulation with nonlinear terms as noted above is used. 


DISCUSSION: 

The desciption of the element is direct. However, it appears that the 
solution procedure is a stepwise linear incremental Lagrangian with a tangent 
stiffness (modified Newton-Eaphson) iteration to avoid cumulative error. The 
total noplinear strain must therefore be available, and the tangent stiffness 
appears to have retained all nonlinearities, ie., through the cubic terms in 
the total displacement. The geometric stiffness matrix does not explicitly 
appear in the equations. 

Slope continuity is imposed at the midpoints of the sides of the element, 
using the Lagrange Multiplier procedure. 

Instability is evaluated as the zero of the determinant of the tangent stiff- 
ness matrix. 

Numerical results are given for static analysis and for instability and 
nonlinear behavior cases, with good accuracy. 
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ELEMENT ID: S-13 

ELEMENT TYPE: SHEBA6 (SHEBA3) 

(General Shell Element 


ASSUMED DISPLACEMENT SHAPE: Complete fifth order 

Eigure ; 


DESCRIPTION: 


Number of nodes : 
Degrees of freedom : 



Yertices 
u u, 




V v,_ 


w w, 


X 


^’y 

XX 

u, 

xy 

u, 

yy 

"’y 

V, 

. XX 

V, 

’xy 

V, 

yy 

w, 

' y 

w, 

XX 

^’xy 

^-yy 

^ nodes, u 

n 


w , n 
n 


REFERENCE: 


1100} 


all 


local direction normal to the edge 
together 63. 

Isotropic or Anisotropic 

Deformation type: stretching, in-plane, shearing, bending 

without shear 

VARIATION OF THIS ELEMENT: 


Eliminations of mid-side-nodes leads to SHEBA3. 


ADVANTAGES OR DISADVANTAGES: 

Advantages : Unachieved accuracy even in case of very coarse grids. 

Successful handling of the most complicated shell problems. 

Eull compatibility. . » 

Disadvantages; Application needs experienced staff. Grid lines in parameter 
plane must be straight. Jumps in thickness provide difficulties. 
Discussion: Based on the concept of natural stress invariants for large rotations 




146 


HANDBOOK OF NONLINEAR FINITE ELEMENTS 


ELEMENT ID: S-14 

ELEMENT TYPE: ARCUS 6 (ARCUS 3) 

(Stiffener layer fully compatible with SHEBA6 (SHEBA 3)) 


ASSUMED DISPLACEMENT SHAPE: Fifth order 



REFERENCE: IlOl] 

VARIATION OF THIS ELEMENT; Free beam is under development (applicable to 
curved grids) 

ADVANTAGES OR DISADVANTAGES: 

Discussion : Based on the concept of natural stress invariants for large 

rotations. 
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ASSUMED DISPLACEMENT SHAPE; Biquadratic Lagrance interpolation for both 
displacements and rotations. 



Derived from, 3D continuum model by "degenerating” assumption in thickness 
direction. 


Number of nodes ; 9 

Degrees of freedom 45 (per node,, i.e., u, v, w, x, 3) 
Deformation type • stretching and in-plane shear; bending 

behavior, stems solely from*--transverse shear. 


REFERENCE; [102] 


VARIATION OF THIS ELEMENT: Possible extension to higher order interpolation 

scheme . 

ADVANTAGES OR DISADVANTAGES; 

Advantages : Good engineering accuracy, simple representation of geometry. 

Disadvantages ; Restriction to analysis of "thick" shells, due to doubtful 
validity of "reduced integration technique" in case of "thin" shells. 

Discussion ; Performal.uce has been investigated within internal research, 

therefore the nonlinear portion of QUABC9 is not incorporated in 

ASKA. 
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ELEMENT ID: C-1 

ELEMENT TYPE; Tetrahedron 


ASSUMED DISPLACEMENT SHAPE: 

6L= (- 1 ^ 

V = (cl~ htj t 2 + 

ClJ= x -^ 



^X;U 


DESCRIPTION: 

This element is an isotropic tetrahedron with 4 nodes and 3 DOF (u,v,w) 
per node. It accounts for material nonlinearities, but does not consider 
geometric nonlinearities. A three-dimensional state of constant stress is obtained 


REFERENCE: [19] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES: 



BEPRODUCIBILITY OF Tim 
ORIGINAL PAGE IS POOP. 


HANDBOOK (CONTINUED) 


ELEMENT ID: C-1 

STRAIN DISPLACEMENT EQUATIONS: 

oCT ^ pf 

e - elastic 
pt - plastic 








z>tpT 


DISCUSSION; 


This derivation considers material non linearities for an isotropic 
tetrahedron. Plastic strains are converted to a column matrix of 
plastic loads. 
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ELEMENT ID: C-2 

ELEMENT TYPE: Triangular, 

Toroidal Con- 
tinuum Element 



ASSUMED DISPLACEMENT SHAPE: 


-A Ll 

/\ -'N 


DESCRIPTION: 

This element is an isotropic toroidal ring element developed for a large dis- 
placement, moderate rotation, elastic-plastic dynamic problem. It has 3 ring 
nodes and 2 DOF (u ,u ) per node. Convected coordinates are used to derive 
the elemental matrices. A two dimensional deformation state is assumed. Both 
geometric and material nonlinearities are included. 

REFERENCE; [54] 


VARIATION OF THIS ELEMENT: 


ADVANTAGES OR DISADVANTAGES; 

The method is stated as being computationally efficient. 
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ELEMENT ID: C-2 


STRAIN DISPLACEMENT EQUATIONS: 


rt ) 


0 


0 


1 , A I 

1 

XrJrj 




is assumed fo be linear. 


DISCUSSION: 

This element is stated as being able to account for materia! non linearities. 

Because convected coordinates are used^ the geometric non linearities that 
arise are accounted for by transformations between the con/ected and global 
systems. The strains are linearly related to deformation displacements relative 
to the convected coordinate system. The dynamic equations of equilibrium are 
solved by a central difference explicit integration procedure which uses total nodal 
forces. 


7,0 REFERENCES 


lo M. J. Turner, E, H. Dill, H. C. Martin, and R. J. Melosh; "Large Deflections 
of Structures Subjected to Heating and External Loads", Journal of the Aerospace 
Scl,, Feb. 1960, p. 97. 

2. H. C. Martin; "On the Derivation of Stiffness Matrices for the Analysis of 
Large Deflection and Stability Problems", Proceedings of Conference on 
Matrix Methods in Structural Mechanics, Dec. 1965, WPAFB, Dayton, Ohio, 

3. B, E. Greene; "Buckling Loads for Columns of Variable Section", SARM 12^ 

The Boeing Company, Seattle, Washington, June 1960. 

4. T. Y. Yang, "Matrix Displacement Solution to Elastica Problems of Beams and 
Frames", Int. J. Solids Str., 1973, Vol, 9, p 829. 

5. T. Y. Yang and R. J. Wagner; "Snap Buckling of Nonlinearly-Elastic Finite 
Element Bars", Computers and Structures, Vol. 3, 1973, pp 1473-1481. 

6. R. T. Haftka, R. H. MallettandW. Nachbar; "A Koi ter Type Method for Finite 
Element Analysis of Nonlinear Struc^Jral Behavior", AFFDL-TR -70-130, Vol. 1 , 
Nov. 1970. 

7. Ro W. H. Wu, and E. A. Witmer; "Finite-Element Analysis of Ijarge Elastic- 
Plastic Transient Deformations of Simple Structures", AlAA Journal, Sept. 1971, 

Vol. 9, No. 9, p 1719. 

8. A. B. Sabir and A. C, Lock; "Large Deflection, Geometrically Nonlinear Finite 
Element Analysis of Circular Arches", Int, J. Mech. Sci., 1973, Vol. 15, p37. 

9. D. G. Ashwell, A. B. Sabir, and T, M. Roberts; "Further Studies in the Appli- 
cation of Curved Finite Elements to Circular Arches", Int, J. Mech, Sci., Vol. 

13, 1971, p 507. 

10. Tc M. Roberts and D, G. Ashwell; "The Use of Finite Element Mid- Increment 
Stiffness Matrices in the Post -Buckling Analysis of Imperfect Structures", 

Int. J. Solids Structures, Vol. 7, 1971, p 805, 

11. R. H, Mallett and P. V. Marcal; "Finite Element Analysis of Nonlinear Structures", 
Proc, ASCE, ST 9, Sept. 1968, p 2081. 

12. R. W. H. Wu, and E. A. Witmer; "Nonlinear Transient Responses of Structures by 
the Spatial Finite -Element Method", A|AA J. Vol. 11, No. 8, Aug. 1973, p 1110. 


153 


13. W. L. Salus, C. Ip, and J. W. Vanderlmden; "Design ConsideraHons of Elastic- 
Plastic Structures Subjected to Dynamic Loads", AlAA/ASME 11th Structures, 
Structural Dynamics, and Materials Conference, April 1970, Denver, Colorado, 

14. T. Belytschko and B. J. Hsieh; "Nonlinear Transient Finite Element Analysts 
With Convected Coordinates", Int. J. for Numerical Meth. in Eng., Vol. 7, 

1963, pp 255-271. 

15. R. H. Gallagher and J. Padlog; "Discrete Element Approach to Structural in- 
stability Analysis", AlAA Journal, June 1963, p 1437. 

16. L. A. Schmit, F. K. Bogner and R. L. Fox; "Finite Deflection Structural Analysis 
Using Plate and Shell Discrete Elements", AlAA Journal, Vol. 6, No. 5, May 

1 968, p 781 . 

17. K. K. Kapur and B. J. Hartz; "Stability of Plates Using the Finite Element Method", 
Proc. ASCE. EM 2, April 1966, p 177. 

18. B. J. Hartz; "Matrix Formulation of Structural Stability Problems", ASCE Jour., 

ST 6, Dec, 1965, p 141. 

19. R. H. Gallagher, J. Padlog and P. P. Bijlaard; "Stress Analysis of Heated Com- 
plex Shapes", ARS Jour., May 1962, p 700. 

20. R. J. Melosh; "A Flat Triangular Shell Element Stiffness Matrix", Proc. Conf. 
on Matrix Methods in Str. Mech., AFFDL-TR-66-80, Dec, 1965, p 503. 

21 . "Derivation of Elements Contained in the ASTRA Structural Program", Boeing 
Document D2-125179-3, Feb. 1972, 

22. J. S. Przemieniecki; "Finite Element Structural Analysis of Local Instability", 

A IAA/ASME/SAE 13th Structures, Structural Dynamics, and Materials Conf., 

San Antonio, Tex., April 1972, TR 66-80; also AlAA J., Vol. 11 , No. 1, 

January 1973. 

23. J. S. Przemieniecki; "Matrix Analysis of Local Instability in Plates Stiffened 
Panels and Columns", Int. J. for Numerical Meth. in Eng., Vol. 5 1972, p 209. 

24. T. Yc Yang; "A Finite Element Procedure for Large Deflection Analysis of Plates 
with Initial Deflections", AlAA J., Vol. 9, No. 8, August 1970, p 1468. 

25. F. K. Bogner, R. L. Fox and L. A. Schmit; "The Generation of Inter-EIement- 
Compatible Stiffness and Mass Matrices by the Use of Interpolation Formulas", 

Proc. Conf. on Matrix Meth. in Str. Mech., AFFDL-TR-66-80, Dec. 1965, 

p 397. 


154 


4 


26. C. Brebbla and J, Connor; "Geometrically Nonlinear Finite Element Analysis", 

Jour. ASCE, EM2, April 1969, p 463. 

27. B. E. Greene; "Stiffness Matrix for Bending of a Rectangular Plate Element 
with Initial Membrane Stresses", SARM 45, The Boeing Company, Seattle, 

Washington, August 1962, 

28. B. W. Papenfuss; "Lateral Plate Deflection by Stiffness Matrix Method with Applica- 
tion to a Marquee", Masters Thesis, C. E, Dept., Univ. of Washington, 1959. 

29. T. Kawai and N. Yoshimura; "Analysis of Large Deflection of Plates by the 
Finite Element Method", Int. Jour, for Numerical Meth. in Eng., Vol, 1, 1969, 
pp 123 -133. 

30. G. C. Pardoen; "Static, Vibration and Buckling Analysis of Axisymmetric Circu- 
lar Plates Using Finite Elements", Computers and Structures, Vol. 3, 1973, 

pp 355-375. 

31 . D. W. Murray and E. L. Wilson; "Finite-Element Postbuckling Analysis of Thin 
Elastic Plates"; AlAA Jour., Vol. 7, No. 10, Oct. 1969, p 1915. 

32. R. W. Clough and J. L. Tocher; "Finite -Element Stiffness Matrices for Analysis 
of Plate Bending"; TR -66-80, 1966, Air Force Flight Dynamics Lab., pp 515-545. 

33. 0„ C. Zienkiewicz; "The Finite Element in Structural and Continuum Mechanics", 
McGraw-Hill, London, 1967. 

34. R. W. Clough and J. L. Tocher; "Finite Element Stiffness Matrices for the Analysis 

of Plate Bending", Proc. 1st Conf. Matrix Method Struc. Mech., AFFDL-TR“66-1 50, 
pp 399-440, 1965. 

35. P. G. Bergan and R. W. Clough; "Large Deflection Analysis of Plates and Shallow 
Shells Using the Finite Element Method", Int. Jour, for Numerical Methods in 
Eng., Vol, 5, pp 543-556, 1973. 

36. R. J. Melosh; "Basis for Derivation of Matrices for the Direct Stiffness Method", 
AlAA, Vol. 1, July 1963, pp 1631-1637. 

37. T. E. Lang and B, J. Hartz; "Finite Element Matrix Formulation of Post-Buckling 
Stability and Imperfection Sensitivity", High Speed Comp, of Elastic Structures, 

Proc. of Symposium, Vol. 2, Liege, Belgium, Aug. 23-28, 1970, p 727. 

38. H. G. Harris and A. B. Pifko; "Elastic-Plastic Buckling of Stiffened Rect. Plates", 
Symposium on Application of Finite Element Methods in Civil Engineering, 

Vanderbilt Univ. , Nov. 13-14, 1969, Ed. W. H. Rowan and R. M. Hackett. 


155 


39. W. G. Carson and R. E. Newton; "Plate Buckling Analysis Using a Fully 
Compatible Finite Element", AlAA Jour., Vol. 7, No. 3, pp 527-529, 

March 1969. 

40. Ao Pifko and G. Isakson; "A Finite-Element Method for the Plastic Buckling 
Analysis of Plates", AlAA Jour,, Vol. 7, No. 1(), Oct. 1969. 

41. R. W. Clough and C. A. Felippa; "A Refined Quadrilateral Element for Analysis 
of Plate Bending", Proc. of 2nd Conf. on Matrix Meth, in Structural Mech., 
AFFDL-tR-68-150, Dec. 1969. 

42. E, P. Popov and S. Yaghmai; "Linear and Nonlinear Static Analysis of Axisym- 
metrically Loaded Thin Shells of Revolution", International Conf. on Pressure 
Vessel Technology, Delft, Netherlands, -Sept-Oct. 1969. 

43. G. P. Bazeley, Y. K. Cheung, B. M. Irons and O. C. Zienkiewicz; "Tri- 
angular Elements in Plate Bending - Conforming and Non -Conforming Solutions", 
Proc. 1st Conf. Matrix Meth. Struct. Mech., AFFD L-TR-66-150, pi547, 1965. 

44. R. G. Anderson, B. M. Irons and O. C. Zienkiewicz; "Vibration and Stability 

of Plates Using Finite Elements", Int. J. Solids Structures, Vol. 4, pp 1031-1055, 

1968. 

45. R. H, Gallagher, R. A. Gellatly, J, Padlog and R. H. Mallett; "A Discrete 
Element Procedure for Thin -She 1 1 Instability Analysis", AlAA Journ. Vol. 5, 

No. 1, Jan., 1967, p 138. 

46. P. V. Marcal; "Large Deflection Analysis of Elastic-Plastic Shells of Revolution", 
TOth Structures, Str. Dynamics, and Materials Conf . , New Orleans, La., April 

1969. 

47. P. V. Marcal; "Large Deflection Analysis of Elastic-Plastic Shells of Revolution", 
AlAA Jour., Vol. 8, No. 9, Sept. 1970, p 1627. 

48. H. Armen, A. Pifko and H. S. Levine; "Finite Element Analysis of Structures in 
the Plastic Range", NASA-CR-1649, Langley Research Center, Feb. 1971. 

49. P. Sharifi and E. P. Popov; "Nonlinear Finite Element Analysis of Sandwich 
Shells of Revolution", Paper 72-356 presented at the 13th Structures, Structural 
Dynamics and Materials Conf., San Antonio, Texas, April 10, 1972, 

50. J. A. Stricklin, W, E. Haisler, H. R. MacDougall and F. J. Stebbins; "Nonlinear 
Analysis of Shells of Revolution by the Matrix Displacement Method", AlAA Jour., 
Vol. 6, No. 12, Dec. ]968, p 2306. 

51 . P. E. Grafton and D. R. Strome; "Direct Stiffness Method Analysis of Shells 
of Revolution Utilizing Curved Elements", AlAA Jour., Vol. 1, No. 10, 

Oct. 1963, pp 2343-2347, 


156 


52. 


J. H. Percy, T. H. H. Plan, S. Klein and D. R, Navaratna; "Application 
of Matrix Displacement Method to Linear Elastic Analysis of Shells of Revolu- 
tion";, AlAA Jour., Vol. 3, No. 11, Nov. 1965, pp 2138-2145. 

53. W. E, Haisler and J. A. Stricklin; "Nonlinear Finite Element Analysis Including 
Higher Order Strain Energy Terms", AlAA Jour., Vol. 8, No. 6, June 1970, 

pp 1158-1159. 

54. T. Belytschko and B. J. Hsieh; "Non linear Transient Analysis of Shells and 
Solids of Revolution by Convected Elements", AIAA/ASME/SAE I4th Structures, 
Structural Dynamics and Materials Conference, Williamsburg, Va., March 1973, 
Paper 73-359. 

55. T. Y. Yang and H. W, Kim; "Vibration and Buckling of Shells Under Initial 
Stresses", AlAA Jour., Vol. 11, No. 11, Nov. 1973, p 1525. 

56. T. Y. Yang; "High Order Rectangular Shallow Shell Finite Element", ASCE 
EMI, Feb. 1973, pp 157-181. 

57. J. L. Sanders; "An Improved First-Approximation Theory for Thin Shells", 

NASA TR R-24, 1959. 

58. R. H. Gallagher, T. Y. Yang; "Elastic Instability Predictions for Doubly Curved 
Shells", Proc. of 2nd Conf. on Matrix Methods in Str. Mech., AFFDL-TR-68-150, 
Dec. 1969. 

59. V. Novozhiiov; "The Theory of Thin Shells", P. Noordhoff Ltd., The Netherlands, 
1959. 

60. D. R. Navaratna, T. H. H. Plan and E. A. Witmer; "Stability Analysis of Shells 
of Revolution by the Finite -Element Method", AlAA Jour., Vol. 6, No. 2, 

Feb. 1968, pp 355-361. 

61. J. A. Stricklin, J. E. Martinez, J. R. Tillerson, J. H. Hong and W. E. Haisler; 
"Nonlinear Dynamic Analysis of Shells of Revolution by Matrix Displacement 
Method", AlAA Jour., Vol. 9, No. 4, April 1971, pp 629-636. 

62. J. Connor and N. Morin; "Perturbation Techniques in the Analysis of Geometrically 
Nonlinear Shells", High Speed Computing of Elastic Structures, Proc. of Symposium, 
Liege, Belgium, Aug. 23-28, 1970, Vol. 2, pp 683-705. 


157 



63. J. H. Argyris "Continua and Discontinua" Proceedings of Conference 
on Matrix Methods in Structu*‘‘:| Mechanics. AFFDL-TR-66-80, Dec, 

1965, pp. 11-190. 

64. A. Pifko, H. S. Levine, H. Armen, Jr., "PLANS - A Finite Element 
Program for Nonlinear Analysis of Structures, Volume I - Theoretical 
Manual" NASA Contract NAS-1-10087, to be published. 

65. W, L. Salus, "A Linear Strain Triangular Element", Boeing Document 
D2-125251-1, The Boeing Co., Seattle, Wash. Dec. 1966. 

66. L. B. Kiersky, "COSMOS, A Computer Program for Structural Analysis", 
Boeing Document D2-4513. The Boeing Co., Seattle, Wash. Jan. 24, 
1962. 

67. H. Levine, at al., "Nonlinear Behavior of Shells of Revolution 
Under Cyclic Loading", J. Computers and Structures, Vol. 3, pp . 
589-617, 1973. 

68. T. Y. Yong and G. Lianir, "Large Displacement Analysis of Visco- 
elastic Beams and Frames by the Finite-Element Method", Jour. 

Appl. Mech. Vol. 41, No, 3, pp. 635-640, Sept. 1974. 

69. K-J. Bathe, H. Ogdemir, E. L. Wilson, "Static and Dynamic Geometric 
and Material Nonlinear Analysis", Report UG SESM 74-4, University 
of California, Berkeley, Calif. Feb. 1974. 

70. R. Sundaramoothy and D. W. Murray, "Completed Local Buckling in 
Wide-Flange Beam Columns," Journal of the Structural Division, 

ASCE, Vol. 99, No. ST6, 1973, Proc. Paper 9774, June 1973, pp . 
1003-1023. 

71. R. Sundaramoothy and D. W. Murray, "Finite Element Solution of 
Inelastic Beam Equations," Journal of the Structural Division, 

ASCE, Vol. 99, No. ST6, Proc. Paper 9773, June 1973, pp. 1025-1041. 

72. M. Epstein and D. W. Murray, "Large Deformation In-Plane Analysis 
of Elastic Beams," to be submitted to Computers and Structures, 
Pergamon Press. Personal correspondance D. W. Murray to R. E. Jones, 
Jan. 2, 1975. 

73. Personal correspondance T. Balderes to R. E. Jones, Dec. 18, 1974. 

74. R. Rung, P. Mason and J. Rosenbaum, "Unsyrametric Analysis of Shells 
of Revolution," Grumman Advanced Development Report ADR 02-11-65.3. 

75. N. Prince, R. Rung, and V. Svalbonas, "Unsymmetric Non-Linear First 
Order Analysis of Layered Orthotropic Shells of Revolution," 

Grumman Advanced Development Report ADR 02-11-66.2. 


158 


76. 




t,sR jE*isa rtiliji'i. . w., — 


t. 


77. 


78. 


79. 


80. 


81. 


82. 


83. 


84. 


85. 


86 . 


87. 


V. Svalbonas, "Unsymmetric Non-Linear Second Order Analysis of, 
Orthotropic Shells of Revolution,” Grumman Advanced Development 
Report ADR 02-11-66.1. 

P. Mason, R. Rung, J. Rosenbaum and R. Ebrus, "Non-Linear 
Numerical Analysis of Axisymmetrically Loaded Arbitrary Shells 
of Revolution,” AIAA Journal, Vol. 3, No. 7, July 1965, pp. 

1307-1313. 

■y. Svalbonas and T. Balderes, "Buckling and Vibration Analysis 
for Stiffened Orthotropic Shells of Revolution,” AIAA Journal, 

Vol. 10, No. 7, July 1972. 

H. S. Levine and V. Svalbonas, "Inelastic, Nonlinear Analysis of 
Stiffened Shells of Revolution by Numerical Integration," ASMS 
Journal of Pressure Vessel Technology, Paper No. 74-PVP-35, May 
1974. 

T, Balderes and A. E. Armenakas, "Free Vibrations of Ring Stiffened 
Toroidal Shells," AIAA Journal, Vol. II, No. 12, Dec. 1973, pp . 
1637-1644. 

P. Sharifi and E. P. Popov, "Nonlinear Buckling Analysis of 
Sandwich Shells," Shell Structures and Climatic Influences. Pro- 
ceedings of Symposium held at the University of Calgary, July 3-6, 
1972. 

P. K. Larsen and E. P. Popov, "Large Displacement Analysis of Visco- 
elastic Shells of Revolution," Computer Methods in Applied Mechanics 
and Engineering 3 (1974) pp. 237-253. 

S. Ahmed, B. M. Irons and 0. C. Rienkiewicz, "Curved Thick Shell 
and Membrane Elements with Particular Reference to Axisymraetric 
Problems," Proc. 2nd Conf . on Matrix Meth. in Structral Mech. 
AFFDL-TR-68-150 (1968) pp. 539-572. 

S. Nagarajan and E. P. Popov, "Nonlinear Elastic-Viscoplastic 
Analysis Using the Finite Element Method," Proc. of Conf. on Com- 
puter Methods in Nonlinear Mech. , Univ. of Texas, Austin, Texas 
Sept. 23, 1974. 

E. P. Popov and P. Sharifi, "A Refined Curved Element for Thin 
Shells of Revolution," Int, Jour, for Numerical Meth. in Eng. 

Vol. 3, 1971, pp. 495-508, 

T. M. Hrudey, "The Application of a Refined Shallow Shell Finite 
Element to Geometrically Nonlinear Structures," Canadian Aeronautics 
and Space Institute Transactions, Vol. 6, No. 2, Sept. 1973, pp. 
92-96, 

Personal correspondance M. J. Jacobson to R. E, Jones, Dec. 16, 1974. 


159 


J. Brass, J. R. Yamame, and M. J. Jacobson, "Effects of Internal 
Blast on Combat Aircraft Structures," Technical Report AFFDL-TR- 
73-136, Vol. Ij Engineers Manual, Jan. 1974. 

Personal correspondance T. Belytschko to R'. E. Jones, Dec. 5, 1974. 

T. Belytchko and A. Marchertas, "Nonlinear Finite Element for 
Plates and Its Application to Dynamic Response of Reactor Fuel 
Subassemblies," J. of Press. Vessel Technology, ASME, Vol. 96, 

Ser. J. No. 4, Nov. 1974, pp. 251-257. 

A.Razzaque, "Finite Element Analysis of Plates and Shells," 

Ph.D. Thesis, University of Wales, 1972 i 

J. H. Argyris, "Recent Advances in Matrix Methods of Structural 
Analysis," Pergamon Press, 1964. 

J. H. Argyris and D. W. Scharpf., "Large Deflection Analysis of 
Prestressed Networks", Journal of the Structural Division, ASCE 
Vol. 98, ST3, March 1972. 

J. H. Argyris, "Three Dimensional Anisotropic and Inhomogeneous 
Elastic Media: Matrix Analysis for Small and Large Displacements", 

Ingenier Archiv, Vol. 34, No. 1, Jan. 1965. 

J. H. Argyris, P. C. Dunne, T, Angelopaubs, B. Bichat, "Large Natural 
Strains and Some Special Difficulties due to Nonlinearity and 
Incompressibility in Finite Elements", Computer Methods in Applied 
Mechanics and Engineering, Vol. 4, No. 2, Sept. 1974., 

J. H. Argyris and D. W. Scharpf, "Some General Considerations on 
the Natural Mode Technique, Part II - Large Displacements", 

T.N. 20, J. Roy. Aer. Soc., Vol. 73, April 1967. 

J. H. Argyris, W. Bosshard, I. Fried, H. M. Hilber, "A Fully 
Compatible Plate Bending Element", ISD Report No. 42, Stuttgart, 
December 1967 . 

H. P. Mlejnek, "Ein Beitrag zur Nichtlinearen Statischen and 
Dynamischen Analyse von Vorgespannten und Rotierenden Platten und 
Schalen mit der Methode der Finiten Elemente", Dr. Ing. Dissertation, 
University of Stuttgart, 1972. ("A Contribution to the Nonlinear 

Static and Dynamic Analysis of Prestressed and Rotating Plates and 
Shells Using the Finite Element Method".) 

J. H. Argyris, K. E. Buch, "A Sequel to Technical Note 14 on the 
TUBA family of Plate Elements", The Aer. Journal of the Roy. Aer. Soc 
Vo. 72, No. 695, Nov. 1968. 


100 . 

101 . 

102 . 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110 V 
111 . 
112 . 


J. H. Argyris, D. W.. Scharpf, ”A Sequel to Technical Note 15: 

The SHEBA Family of Shell Elements for the Matrix Displacement 
Method, Part III - Large Displacements”, J. Roy. Aer. Soc., 

Vol. 72, No. 694, Oct. 1968. 

J. H. Argyris, M. Hoose, G. Malej annakis , etc., to be published. 

G. Maseczek, H. Parisch, "Finite Element Study of Geometrical 
Nonlinearities”, submitted for publication. 

J. A. Stricklin, and N. E. Haisler, "Survey of Solution Procedures 
for Nonlinear Static and Dynamic Analyses". 

K. J. Bathe, E. Ramm, and E. L. Wilson, "Finite Element Formulations 
for Large Deformation Dynamic Analysis", Int. Journal for Numerical 
Methods in Engineering (date not known) . 

K. J. Bathe, E. Raimn, and E. L. Wilson, "Finite Element Formulations 
for Large Displacement and Large Strain Analysis", University of 
California Report UCSESM 73-14, Sept. 1973. 

0. C. Zienkiewicz, and G. C. Nayak, "A General Approach to Problems 
of Large Deformation and Plasticity Using Iso-Parametric Elements", 
Presented at Third Conference on Matrix Methods in Structural 
Mechanics, WPAFB, Oct. 1971. 

J. T. Oden, "Finite Elements of Nonlinear Continue", McGraw Hill, 
1972. 

R. G. Vos, "Development of Solution Techniques for Nonlinear 
Structural Analysis", Final Report NASA Contract NAS8-29625 
Marshall Space Flight Center, Alabama, Sept, 30, 1974. 

S. W. Key, "A Finite Element Procedure for the Large Deformation 
Dynamic Response of Axisymmetric Solids", Computer Methods in 
Applied Mechanics and Engineering, Vol. 4, pp. 195-218, (1974). 

J. T, Oden, and L. C, Wellford, "A Finite Analysis of Shock and 
Vibration Problems - Recent Developments, Shock and Vibration Digest, 
(to appear) . 

L. E, Wellford, and J. T. Oden, "A Theory of Discontinuous Finite 
Element Galerkin Approximations of Shock Waves in Nonlinear Elastic 
Solids", (in review). 

L. C. Wellford, and J. T. Oden, "On Some Finite Element Methods 
for Certain Nonlinear Second-Order Hyperbolic Equations, Sept. 1974 
(in review) . 


161 



113. 

114. 

115. 

116. 

117. 

118. 

119. 

120 . 


J. T. Oden, "Approximations and Numerical Analysis of Finite 
Deformations of Elastic Solids", Nonlinear Elasticity , ed. by 
W. Dickey, Academic Press, New York, 1974, pp. 175-228, 

J, T, Oden, "Formulation and Application of Certain Primal and 
Mixed Finite Element Models of Finite Deformations of Elastic 
Bodies", Lecture Notes in Computer Science, (edited by G, Goos 
and J. Harminas) , Vol, 10, Computer Methods in Applied Sciences 
and Engineering, Part 1, Springer-Verlag, Berlin-Heidelberg, 
pp. 334-365, 1974. 

J. T. Oden, and J. E. Key, "Analysis of Finite Deformations of 
Elastic Solids of the Finite Element Method, Proc, lUTAM 
Colloq. High Speed Computing Elastic Strut. , Liege, 1971. 

J. T. Oden, and J. E. Key, "Numerical Analysis of Finite 
Axisymmetric Deformations of Incompressible Elastic Solids of 
Revolution", International Journal of Solids and Structures, 

Vol. 6, pp. 496-518, 1970. 

J. T. Oden, and T, Sato, "Finite Strains and Displacements of 
Elastic Membranes by the Finite Element Method", International 
Journal of Solids and Structures , Vol. 3, July 1967, pp. 471-488. 

J. T. Oden, and W. K. Kubitza, "Numerical Analysis of Nonlinear 
Pneumatic Structures", Proceedings, International Colloquium on 
Pneumatic Structures , Stuttgart, Germany, pp. 87-107, 1967. 

R. H. Gallagher, and G. R. Thomas, "The Finite Element Mechanics 
in Plate and Shell Justability Analysis',' paper presented at 
Fourth Australian Conference on the Mechanisms of Structures and 
Materials, Brisbane, Australia, Aug. 1973. 

Jerrell M. Thomas "A Finite Element Approach to the Structural 
Instability of Beam Columns, Frames, and Arches," NASA TN D-5782, 
May 1970. 


162 



